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Abstract 

^ I We report here the status of different gauge conditions in the canonical formulation of quantum 

i electrodynamics on light-front surfaces. We start with the massive vector fields as pedagogical 

models where all basic concepts and possible problems manifestly appear. Several gauge choices 
are considered for both the infinite and the finite volume formulation of massless gauge field 
electrodynamics. We obtain the perturbative Feynman rules in the first approach and the quantum 
^ . Hamiltonian for all sectors in the second approach. Different space-time dimensions are discussed 

in all models where they crucially change the physical meaning. Generally, fermions are considered 
as the charged matter fields but also one simple H-1 dimensional model is discussed for scalar fields. 
Finally the perspectives for further research projects are discussed. 



X 



*This paper is a part of a qualifying thesis for habilitation in the Faculty of Physics of the Warsaw University. 



1 



Contents 



I introduction 



II Massive Electrodynamics 



1 Theory without gauge fixing condition 



1.1 Maxwell theory with a mass term 



1.1.1 Vector fields with external currents 



6 
6 

6 

6 

9 

11 

13 

1 Vector field scctorl 14 



1.2 Gross- Treiman model 



1.2.1 Interaction with fermion fields 



2 Lorentz covariant gauge 



2 Interactions with fermions 16 



p. 3 Perturbation tlieory 19 



2.3.1 Field operators in the interaction representation 19 



2.3.2 LF perturbative calculations of the S-matrix elements 22 



3 LF Weyl gauge 



23 

3.1 Vector fields sector 23 



3.2 Free quantum fields 24 



III Light Front QED 



27 



4 Class of LF Weyl gauges 



27 

4.1 Gauge field sector in 1+1 dimensions 27 

4.2 Gauge field sector in D+1 dimensions 29 



5 General axial gauge 



31 

5.1 Gauge fields in 1+1 dimensions 31 



5.2 Gauge fields in 3+1 dimensions 33 



6 Flow covariant gauge 



35 

3.1 Model in 1+1 dimensions 35 

|6.1.1 Gauge field sector 36 

36 



3.1.2 Interaction with fermion fields 



.1.3 Perturbation theory 37 



3.2 Highcr-dimcnsional model 38 



6.2.1 Perturbation theory 41 



Electrodynamics of charged scalar fieldsj 43 
[7.1 LF We~ gauge in 1+1 dimensions] . . 43 



.2 Dressed scalar fields 45 



.3 Perturbation theory] 47 



2 



iV Finite volume QED 



8 LF Weyl gauge 

p.l The DLCQ methodi 

g.2 Canonical formalism for LF Weyl gauge 

^■2.1 Proper zero mode sector . . . . 

^■2.2 Normal mode sector 

^.2.3 Fermion sector 

g.3 Quantum theory 

B.3.1 Translation Generators 

3.3.2 Implementation of Gauss' Law . 



9 Covariant gauge 

p.l Normal mode sector 

p. 2 Proper zero mode sector 

p. 3 Fermion field sector 

p. 4 Quantum theory and physical states 



10 Transverse Coulomb gaug^ 

10.1 Proper zero mode sectoi .... 

10.2 Normal mode sector 

10.3 Fermion field sector 

10.4 Canonical quantization 

10.5 Physical gauge conditions . . . 

V Conclusions and perspectives 



Appendices 



A LF Notation 

A.l Coordinates . 
A. 2 Dirac matrices 



B Green Functions 

B.l Feynman propagator functions! . 
B.2 Integral operators in x~ direction 
B.3 Inverse Laplace operators .... 
B.4 Finite volume Green functions . . 



References 



Noninmensum gloriantes in alienis laboribus 

2 Cor. 10.15 



Part I 

Introduction 

Half a century ago Dirac Q has proposed 3 different forms of relativistic dynamics depending on 
the types of surfaces where independent modes were initiated. The first possibihty when a space-hke 
surface is chosen (named by Dirac instant form) has been used most frequently so far and is usually 
called equal-time quantization (ET). The second choice is to take a surface of a single light wave 
(called by Dirac front form). This kind of quantization will be discussed in this work and is commonly 
referred to as light-front formalism (LF). The third possibility is to take a branch of hyperbolic 
surface x^^x^ = ■, xq > (named by Dirac •point form). Though a single attempt in this last 
approach can hardly be found, we quote also this choice because very frequently the LF formalism is 
erroneously called light-cone quantization - a name which correctly should be connected with a special 
case (when — > 0) of the point form of dynamics Q. Then it took almost 20 years before Dirac's 
idea of front form of dynamics was applied by physicists. 

At an ET surface, any two different points are spatially separated, therefore fields at these points are 
naturally independent quantities. For a LF surface things are different because for any point, there is a 
null-direction and all points lying along this direction are light-like separated, thus fields at these points 
may have nonvanishing commutators. For all other directions on a LF surface, points are spatially 
separated and fields at these points are independent. This property of the LF formalism has been 
used to infer the behaviour of light-cone commutators |^ from the respective LF commutators. 
When Weinberg considered the scalar field theory in the infinite momentum frame (IMF) |^, 
he has found great simplifications of the " old-fashioned" time-ordered perturbation theory. Then 
the close relation between both the IMF formalism @, §, §, g, §, @, p| and the LF 

formalism IQ, |14], [^], [|l^], [18| was established, and these two names were quite frequently 



interchangeably used. Since then, a lot of successful attempts to phenomenology has been done |19] 
most of them being based on the LF perturbative calculations [20|. 



Another possibility started over 10 years ago when Brodsky and Pauli [^] have introduced so-called 
discrete light-cone quantisation (DLCQ), where periodic boundary conditions have been imposed 
on fields in a finite volume of the LF surface. This approach aimed at the nonperturbative approach 
opened a broad scope of both the theoretical and numerical studies. Another important attempts 



in the nonperturbative directions are the renormalization of Hamiltonians [22| and the relativistic 
bound-state problem [^]. Though, at the time, they are beyond the scope of our paper, the next 
steps that one can take starting from our analysis may be made precisely in these directions. 
Main topic of this paper is the canonical description of quantum electrodynamics (QED) within 
the LF approach. This model is best known in the standard ET formulation but its LF analysis 



has also been studied in many papers since the first formulations in [0], [^, [14|. However, one 
important question of gauge fixing conditions, which can be imposed of gauge field potentials, has 
been not solved so far. Nearly all papers use, seemingly the only possible, the gauge fixing condition^] 
= A- = named the Ught-cone gauge (LC-gauge) or the IMF-gauge. Evidently such a choice 
simplifies enormously many problems, specially those which appear in the sectors of charged matter 
fields. Also this gauge allows for the propagation of only physical photons in the sector of gauge fields. 
However, there is a minor inconsistency in the LC-gauge, because the Feynman perturbative rules 



^Notation is given in Appendix A.l 
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found in the ET quantization pi] have the causal Mandelstam-Leibbrandt (ML) p^] prescription 



(n ■ k + iesgn{n* ■ k)Y 



for the spurious poles of gauge field propagator, while the corresponding rules found in the LF quan- 
tization give the Cauchy Principal Value (CPV) prescription for these poles. While for Abelian 
models both prescriptions are equally correct then non-Abelian models need exclusively the ML- 



prescription [26|. This important caution should be kept in mind during the present discussion where 
Abelian interactions don't falsify the CPV prescription but only prepare a future developments with 
non-Abelian interactions. 

In the canonical LF approach one chooses one null-direction, usually parameter of dynamical 

evolution and treats all 5+ operators as 'time' derivatives. Therefore the ET and the LF descrip- 
tions of the same relativistic models formally may look quite different, specially the structure of their 
canonical conjugate momenta. Here we will not dogmatically follow Dirac's general method of canon- 
ical description for constraint system |27|. But to the contrary, we will rather focus our attention on 
important physical points. First, we will consider as truly constrained variables only those fields which 
satisfy non-dynamical equations (i.e. with no 9+ terms). Such equations can be formally solved for 
non-dynamical fields and therefore, non-dynamical and dependent variables are understood as syn- 
onyms here. By definition, all other fields are considered as dynamical and independent variables. 
Second, canonical momenta for dynamical fields are introduced and canonical Hamiltonians are de- 
fined via the Legendre transformation. If necessary, the terms which contain 9+ will be expressed by 
means of corresponding canonical momenta and only in this way the conjugate momenta can appear in 
the canonical LF Hamiltonian. For all variables which appear in the Hamiltonian we have dynamical 
equations of motion which are of the first order with respect to 5+, and we demand that these Hamil- 
ton equations follow from the canonical Hamiltonian by means of the classical Poisson-Dirac brackets. 
Next, having a consistent classical canonical structure one can define a consistent quantum theory 
where Poisson-Dirac brackets are replaced by quantum commutation relations, and the Heisenberg 
equations for quantum field operators have the same form (up to some ordering of noncommuting 
operators) as the corresponding classical counterparts. Third, for the interacting system of gauge and 
matter fields one can carry the above canonical method in steps, starting with the gauge field sector 
where all interactions with matter are described by linear couplings to arbitrary external currents. 
Having eliminated all gauge constraints one can write an equivalent Lagrangian which describes the 
gauge sector by means of fewer unconstrained fields. Finally the canonical quantization of complete 
system can be carried out quite easily having the constraints for gauge and matter fields transparently 
separated and solved.^ 

Our paper is organized as follows. In part II a pedagogical model of massive vector electrodynamics 
with fermion matter fields is analysed. Though this is not a true gauge field model, it is discussed 
in three cases of no gauge fixing, LF Weyl and Lorentz gauges. In part III the QED with fermions 
is presented for several different gauges, and perturbation Feynman rules for practical S-matrix cal- 
culations are given in each case. Also QED for charged scalar fields is presented in 1-1-1 dimensions 
for the LF Weyl gauge. In part IV we analyse the DLCQ approach to QED with periodic boundary 
conditions for gauge fields and antiperiodic conditions for fermion fields imposed on the finite volume 
LF. Also here several gauge conditions are discussed and the canonical analysis is carried out subse- 
quently in different sectors for various gauge field modes and fermion fields. In Conclusions all results 
are generally discussed special attention being paid to the further developments. The Appendices 
contain all notations and definitions of different functions used throughout the main text. 
Most of these results have not been published so far, otherwise due references to original papers are 
given in the paper. 

■^This method of quantization is similar to that proposed by Jackiw [pq]. 
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Part II 

Massive Electrodynamics 



1 Theory without gauge fixing condition 

Our first model, which we analyse along the lines indicated in the Introduction, will be the electro- 
dynamics of massive vector field -B^ with the mass term m^. Its kinetic term has the Maxwell form 
built with Bf^i, = d^B^ — di,B^ and it couples with fermion fields via the electrodynamic current 
term gip'j^'ipB^. Though this model has been consistently quantized on LF by Yan over 25 years ago 
via the Schwinger action principle |^9|, we decided to derive his results in our original method 
of quantization for constrained systems. This model is even quite a pedagogical example because it 
contains all obstacles and restrictions which later will be encountered in more physical cases. 

1.1 Maxwell theory with a mass term 

We start our analysis with the following Lagrangian density: 

1 TnP — 

{d+Bi - diB+){d.Bi - diB_) + \{d+B_ - d+B_f - -^{diBj - djBif 

2 

^-{2B+B^ - Bf) + V;(i5^7^ - eB^^^' + M)^ , (1.1) 

where the explicit LF coordinates are introduced for the vector fields. It is known that only a half 
of the fermion degrees of freedom are dynamical fields, therefore in order to separate the vector field 
constraints from the fermion field ones, we have decided to start with a subsystem where the fermion 
contribution is solely described by external currents j'^ coupled linearly with B^. Because in such 
a simplified model the dynamics of fermions is omitted, one cannot argue in favour of the external 
current conservation law d^j^ = 0. Even more we stress that all components of the current should be 
taken as independent classical quantities which allow for taking functional derivatives with respect to 
them^ 

For such a simplified model one can easier find a set of dynamical independent field variables which 
identically satisfy the constraints and have correct dynamical properties. Later we can write an equiv- 
alent Lagrangian with these independent modes and external currents j^, and then by reintroducing 
fermion kinetic terms we end up with the complete equivalent Lagrangian and finally, quantize a 
complete system without any effort. However this scenario is valid only if the condition of mutual 
independence of fermion and vector fields is satisfied. Otherwise some obstacles will appear and (for- 
tunately or not) this situation will appear in our first model. 

1.1.1 Vector fields with external currents 

Our first task will be to analyse a subsystem of vector fields, so we start with the Lagrangian density 
= (d+Bi - diB+){d.Bi - diB^) + \{d+B^ - d^B+f - ^{OiBj - d.Bif 

^This property of external currents is crucial in checking the consistency of the simplified model with the one we 
have started with. 



+ 
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m 



+ ■^{2B+B^-Bf) + B^f , (1.3) 

and taking the variations of vector fields B^ we generate tlie Euler-Lagrange equations 

d+{d+B^ -d^B+- diBi) = {jr? - A^)B+ + j- , (1.4) 

-a_(a+5_ - + diB,) = (m^ - a^)s_ + j+ , (1.5) 

{2d+d- -A±+ rr?)Bi = d^{d+B^ + d-B+ - djBj) + f. (1.6) 
As tlie consistency condition for these equations we have 



d+B^ + d-B+ - d,B, = — ^ + + dd 



(1.7) 



so, together with (|l.4| ), we end up with two independent equations with d^B^ - but the excess of 
dynamical equations indicates the presence of constraints. This observation will become transparent 
after introducing the canonical momentum IT". Generally at LF, one introduces canonical momenta 
n'' as the functional derivatives of Lagrangian by d+B^, respectively. In our present case this gives 



n 

n+ 



d+B- - d-B^ 
d-Bi - diB_ 
, 



(1.8) 
(1.9) 
(1.10) 



which according to Dirac's nomenclature indicates the presence of primary constraints (|1.9| ) and ( |1.10| ) 
because they do not determine d+Bi and c?+i?+. Constraint 11+ is characteristic also for the gauge 
fields and here the component B^ is a dependent variable which can be determined from Eq. (|l.7] ) as0 



B: 



2d. 



n- - diBi + ^2 ( d+r + d-j- + dif 



(1.11) 



Other constraints 11* — d-Bi + diB- ~ regularly appear on LF for covariant relativistic fields and 
they can be ignored only after the canonical Hamiltonian is calculated 



njcan 

rtn 



U'd+B^ +U'd+Bi- Cb 

i(n-)2 + ^{diBj - djBif + ^m^Bf - B^j- - Bif 



(1.12) 



+ B, 



~d-{Tr + diBi) + (Ax - m^)B. 



J 



We notice that the expression in the square bracket identically vanishes due to Eq.(L5) which deter- 
mines the field component i?_ 



5_ 



1 



[d^{ir + d^B,) + 3+] . 



A_|_ — 

For other dynamical field components we have the equations of motion 



(2d+d- - Ax + m^) n- = -d+i+ + d-j 

A 

-.2 



(2d+d- 



A, S 



d+r + d-j- + d^f] +j 



(1.13) 

(1.14) 
(1.15) 



^Integral operators and their convolutions which appear hereafter are introduced in Appendix 



B.3 
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which contain the troublesome expression Here, contrary to the complete theory with fermion 

dynamics, we cannot impose the external current conservation for curing this problem but rather 
introduce new field variables 



n = + 

Bi = Bi + ■ 
m 

which satisfy the dynamical equations 

(id+d- - + m^) n = d-j- + (Aj 

(29+ a 



1 



2a 

m? 2d- 



m 



2d- 



L - A I + B, 



m 



2d- 



■3 



+ f 



(1.16) 
(1.17) 

(1.18) 
(1.19) 



Now the independent dynamical equations have the form of Hamilton equations and the Hamiltonian 
density can be rewritten in terms of these new fields 



l(n 

2^ 



1 



2d- 
1 



*d^{u + diBi 



djBif + \? 



Bi 



2m^ 



B, 



_ AJ_ 

2d- 

d, 1 



m? 2d- 



* J 



(1.20) 



Next, one can demand that Eqs.( |1.18| and |1.19| ) follow canonically from this Hamiltonian and this 
produces Dirac brackets^ at LF 



2dPL {n(x+,f),n(x+,y)} 
2d'L {n(j;+,x),Si(x+,y)} 

2d^_[Bi{x+,x),Bj{x^,y)] 



-(A_L - m^)5'^{x-y) 



DB 
DB 

DB 



, 



djdj 



(1.21) 
(1.22) 

(1.23) 



In this way we have found a canonical structure with Hamiltonian and brackets which is almost 
ready for further canonical quantization, however, before we will make this next step, we should 
check its consistence with the primary Lagrangian ( |1.3D . When we functionally differentiate the 
primary Lagrangian with respect to the external currents j^, we obtain primary fields S^. Similarly, 
differentiation of the consistent canonical Hamiltonian should produce —B^ but in this case Eq.( 1.2d| ) 
do not correctly reproduce the expression for B+ (pip . Physically this means that infiuence of the 
fermion fields on the massive vector fields cannot be reduced here to the presence of charged matter 
currents, but contrary, these two kinds of fields have non-vanishing mixed brackets at LF. Keeping 
this limitation in mind one can yet introduce the effective Lagrangian for vector field sector. First we 
parameterize the primary vector fields as 



B+ 
B_ 
Bi 



d+ h 7^ * 



d. 



m 



m 



d. 



mcj) + diCi + ^ * j" 
d- 



di- + Ci , 
m 



(1.24) 

(1.25) 
(1.26) 



^Though in our analysis we didn't follow Dirac's procedure, we have decided to call these brackets Dirac ones in 
order to stress their consistency with all constraints present in the primary description of the model. 
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where new fields (j) and Cj satisfy equations of motion 

(29+a_ - + m^) ,^ = -1 + - ddi) - * j+ , (1.27) 

2a+a_ - Ax + m^) a = / + (1.28) 

When this parameterization is introduced into the primary Lagrangian (p..3|), we end up with the 
expression 

+ d+(t)d-(t) - - [dict)) - —4) - o K- * J + 3 + J + — J 

2 2 2 \0- J m m m 

+ Cif - (mcl) + d^Ci)j- * j+ , (1.29) 

where evidently fields Ci and 4> are independent modes. We stress that now there are no constraints 
in this effective Lagrangian and one can reintroduce the fermion kinetic terms and substitute external 



currents by the fermion currents. Next the canonical quantization leads to the same results as in |11], 
especially the mixed commutators for the (p field and the fermion fields are nonzero. 

1.2 Gross-Treiman model 

More than 25 years ago, Gross and Treiman have proposed a modified model for describing the 
interaction of massive vector fields with fermions ||3^. They have been specially interested in the 
commutators with a smoother behaviour near the light-cone, and their model can be described by 
the following modified Lagrangian: 

Cgt = {d+v, - a,y+)(a_Fi - diV.) + i(a+y_ - d+v^f - ]{diV, - d^Vif (i.so) 



2v ^ - ^ 4 



where another scalar field has been added with apparently incorrect signs in quadratic terms. Now 
we can take the sector of vector and scalar fields by omitting the fermion dynamics 

4^ = {d^Yi - diY^)(d^Y, - diY^) + i(a+y_ - d^Y^f - \(diY^ - d^Yif 
and then we change variables 

y, = i?^ + ^. (1.32) 

m 



We notice that the Lagrangian ( |1.31[) splits into two parts £j^j, = £^+^0) where Lb is the previously 



analysed Lagrangian ( |1.3| ) and describes solely negative metric scalar field 

£^ = + ^'A' + (1-33) 
2 2 m 

Therefore we need to analyze only this last contribution of ^ fields here and then adopt previous 
results for fields. The Euler-Lagrange equations 

(25+5_ - Ax + m2)(/> = -(9+i+ + d^r + ^^f). (1-34) 

m 
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suggest the change of field variables 



2m 9„ 



which removes term from the equation of motion 



m 



2d. 



If we introduce the canonical momentum for the d) field 



m 



(1.35) 



(1.36) 



(1.37) 



then the canonical Hamiltonian density is 

n/can 



di 1 



+ 



1 



m 2(?_ 
1 1 



1 1 

m 2(?_- 



* J 



m \ m2d- 
and we find the Dirac bracket at LF 



H +dij 



2d^_ {^{x+,x), Hx+y)}^^ = S^x - y). 



(1.38) 
(1.39) 



At last we may add the above results to those found previously and we write down the density of 
effective Hamiltonian as 



n 



GT 



2 



-B 



J-d- 



A I — 



m 



* (n + diBi) + £ 

.4. 1 

-J 



2a_ 



A 1 — m? 



m 



n + 9,5, 



(1.40) 



and nonzero Dirac brackets at LF for all independent field variables 

2a^{n(x+,x),n(x+,y)}^^ = (a^- 



2d'-t[B,{x+,x),Bj{x+J)] 
2dl[4>{x),m} 



DB 



DB 



Sij - 
5{x-y). 



^)5\x-y) , 


(1.41) 


y), 


(1.42) 




(1.43) 



Next, we can easily check the consistency of this effective Hamiltonian with the starting point in 
( 1.31 ) - in both cases currents couple linearly to fields^ 



B, + di^ = Bi + di^ , 
m 



m 
m 



1 



B+ + d, 



1 



A_|_ — m? 
A I — 



* (H + diB,) + ^ 
m 



m 2d- 



H + a-R; 



m 



(1.44) 
(1.45) 

(1.46) 



In the case of V+ one needs to use the equation of motion (1.36) in order to have the term dj, 
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From these results we expect that the above structure of brackets wih survive also in the interacting 
theory when complete dynamics of fermions will be reintroduced. In order to achieve this aim most 
easily, we build here the equivalent Lagrangian 



V4> 

GT 



1 



^ * a_n + d+Bi 6ij - didj - 

A I — V A I — m"' 



1 



(1.47) 



which is nonlocal at LF but directly leads to Dirac brackets (1.41, 1.42| , 1.43) and equations of motion 
( pTlSl , |L19| , |L38|) are generated as the Euler-Lagrange equations. 



1.2.1 Interaction with fermion fields 

Having found the effective equivalent prescription of vector and scalar fields, we can add fermion 
kinetic terms and study the complete Gross- Treiman model 

1 „„ 1 



GT 



d^u-— — - * a_n + d+Bi 6ij - didj- ^* ) d^B, - 

A_L — y J J _ ^2 



+ 



+ 



2 ~n 

H 4) 

2 ^ 2 



2 * (n + diBi) + ^ 

A I — m 



m 



Bi + di—A -e\/2V'U4 



29_ 



A I — 



n + diBi 



m 



(1.48) 

where we have introduced two components of fermion fields 'ip± = A±ip^. The components and 
V'l of fermion fields satisfy the non-dynamical Euler-Lagrange equations 



V2[id- -e(a_$)]V'_ 
V2[-^^--e{^-^)]^pl_ 



1 



M7° - idia' + ea' Bi + 



di 



m 



+ idia' + Bi + 



* (n + diB,) + ^ 



a,; 



ra 



m 



(1.49) 
(1.50) 

(1.51) 



where additional notations ^, and 11 are introduced for convenience and clarity. Next we can write 
the formal solutions for these non-dynamical fermion components. 



1 g-ie* 1 ^ l^g^e*^^ 



1 



id- 



1 



id- 



le'^ 



(1.52) 
(1.53) 



In this way there are only dynamical fields left and their quantization is straightforward. The Hamil- 
tonian density has the following form: 



n 



eff 
GT 



2 



2 " 2 



Bi-djBif + '^B, 



1 



L*(e-*^)+eV2V^t_V'+^ 



m 



n + diBi 



m 



(1.54) 



^See Appendix |A.2| for the complete LF notation of the Dirac matrices and the projection operators A± . 
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while the non-vanishing (anti) commutators at LF have the expected form 
29^ [U{x+,x),U{x+,y)] = i{A± - m'^)5^{x~y) , 



Bi{x'^,x),Bj{x~^,y) 



-i 



rm? 



2dl = i5{x-y) , 

{tp\.{x^,x),tp+{x^,y)} = -^A+d{x-y). 



(1.55) 
(1.56) 
(1.57) 
(1.58) 



The quantum theory defined above is formulated in a natural way in the Heisenberg representation, 
where the whole dynamics of the system is connected with the quantum field operators. For pertur- 
bative calculations the interaction picture is a more convenient choice, but here we will not proceed 
in this direction, leaving this important issue to a more physically relevant models which will be 
discussed later. Rather we end up our discussion of the present model with pointing out its two 
interesting properties. The first one is connected with the scalar field which evidently interacts 
with fermion fields 



— (Aj 

m 



m 



ey/2 
'2m 



id- V " J id. 



However, one can define another field 



(1.59) 



(1.60) 



which, due to equations of the fermion fields 



id+ip^ 



1 



2d. 



A, 



m 



m 



1 

+ 2 



- n + diB, 
d 



-id+tp 



1 



2d. 



A, 



m 



n + diBi 



m 
1 



will satisfy the free field equation of motion 

{2d+d- - A^ + m^)^ = 0. 



M7° + idia' + ea' [Bi + ^ 'c, 

m 



(1.62) 



(1.63) 



However, the price for such a free evolution has to be paid at the level of LF commutators where now 
the non-vanishing relations with (j) are 



2dl [<t){x+,x),^{x+,y)\ 
2&t [0(x+,x),V+(x+,y)] 

2d^_ 



i5^{x-y) , 



m 



-il'+{x~^ ,x)5'^{x — y) 



,x),i)\{x'^ ,y)\ = -—ip!^{x^,x)6^{x-y). 



m 



(1.64) 
(1.65) 

(1.66) 
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The second property is connected with the components of the vector V^: 

Vi = Bi + d^ 



m 
1 



A 



* (n + diBi) + 



m 



1 



+ 



2a_ 



A 



m 



m 



- n + diBi 



which satisfy equations of motion 

{23+3- - A^ + m^)V+ = eV2V+V'+ , 
{23+3- - Ax + m2)y_ 



{23+3- - Ax + m^)Vi 



e 

7i 



1 1 

id- id- 
1 

i3- 



and have the only non-vanishing commutation relations 

25f [V^{x+,x),V,{x+,y)] = -y). 



(1.67) 
(1.68) 

(1.69) 

(1.70) 
(1.71) 

(1.72) 

(1.73) 



Therefore can be interpreted as the massive vector field in the so-called Feynman gauge. We see 

that modification of the massive theory, which solves the problems connected with singular behaviour 
of fields on LF, can be interpreted in terms of gauge condition, though for non-vanishing mass ^ 0, 
this model is not a true gauge theory. Therefore we suggest that also other modifications by means 
of truly gauge fixing terms may be worth to be studied in the canonical LF formulation of massive 
electrodynamics . 



2 Lorentz covariant gauge 

In the previous section we have learned that addition of the extra scalar degrees of freedom can change 
LF commutators to the form proper for the Feynman gauge, which is a special case of the Lorentz 
covariant gauges. General covariant gauge is implemented into Lagrangian by means of the scalar 
Lagrange multiplier field A and it takes the form d^B^ = — aA. For a = it becomes the Lorentz 
gauge and for a = 1 it becomes the Feynman gauge. 

As previously, we start with the classical fields theory defined by the Lagrangian density 

CZss = {d+Bi-diB+){d-B,-diB-) + ^{d+B--d-B+f-^{diBj-djBif 
+ (b-B+ - ^Bf^ + {i^^^d^ - e^B^ -M)^ 

+ A{d+B-+d-B+-d^Bi) + ^A^. (2.1) 

We expect that the canonical structure of vector fields i?^, is independent of the fermion fields, therefore 
we will study the sector of vector fields with external arbitrary currents first. 
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2.1 Vector field sector 

Omitting the fermion kinetic terms in and inserting j^^ in the place of etpjil^ we obtain 
^Tri^ass = {d+B,-d^B+){d^B,-d,B.) + ^{d+B_-d.B+f-^{diB,-d,B,f 

+ m^(^B_B+-^Bf^+B^f + A{d+B^ + d-B+-diB^) + ^A\ (2.2) 

and next we generate the Euler-Lagrange equations for all field variables 

(2a+a_ - Ax + m^) 5_ = (l-a)a_A-j+, (2.3) 

{2d+d- - Aj_ + m^) B+ = (1 - a)d+A - j' , (2.4) 

(2a+a_ - Ax + m2) = (1 - a)aiA + / , (2.5) 
d+B^ = -d-B+ + diBi-ak. (2.6) 

The consistency condition for these equations has the form of a dynamical equation for A 

(2d+d- - Ax + am^) A = + d-j' + dif , (2.7) 

and we see that in order to have no imaginary mass tachyon, we must keep the parameter a > 0. 
Just like in the previous section, we have two different equations with (^]^) and ( |2.6D so there 

is a constraint 

2d^ {djBj -aA- d-B+) = (Ax - m^)B_ + (1 - a)d-A - j+. (2.8) 
If we define the canonical momentum conjugated to the field B- 

n- = d+B^ -d-B+ + A (2.9) 

then, using the gauge condition ( |2.6| ) we obtain another constraint 

2d-B+ = diBi - + (1 - a)A. (2.10) 



In the dynamical equation for A (2/7) there is the term and in order to remove it, we redefine 

the Lagrange multiplier field 

\ = A--^*3+. (2.11) 

In this manner we have selected the set of fields (n~, A) which satisfy dynamical equations of 
motion 

(2a+5_ - Ax + am2)A = (Ax - «"^')^ * J+ + d-j~ + ^^f > (2-12) 



(25+5_ - Ax + m^)n- = m\l-a)^X + — *j+j+2^-r+^^3\ (2.13) 
(2a+a_ - Ax + m2)5, = (l-a)5A + / + ^^^*9.i+ , (2.14) 
so they are independent canonical variables. Also we have dependent fields B^ and By. 



B_ = ^ 2 * (n~ + ^« A - 2A) , (2.15) 

A I — 



2d. 



diB, - n- + (1 - a)\ + * r 
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(2.16) 



Though the other canonical momenta apparently represent primary constraints (according to Dirac's 
nomenclature) 

U\x)-d-Ai{x) + d^A_{x)^0 , (2.17) 
n+(x)psO, (2.18) 
Uxix) « , (2.19) 

we notice that they are absent form the canonical Hamiltonian density^ 

Hcan = {d+B^)U- + (5+S,)ff - ^ = ^ (n" - ^ - ^ * i+) ' + ^ (diBj - djBif + "^Bf 

+ + ^ * j^) diBi - j- ^-^—^ * d. {U- + diBi - 2A) - Bif (2.20) 



m 



a 



2d^_ {Bi{x) 


'n-(y)}^s 


= df6^{x- 


y), 


2dl {n-(x) 




= A^5^{x- 


-y) 


2dl {Bi{x 




= -5ij6^{x 


-y) 


2d^ {A(f 




= -df6^{x 


-y) 


2dt{X{x),Xm^^ 


= w?6^{x- 


-y), 


2df_ {\{x) 


'n-(y)}^B 


= rr?6^{x- 


-y), 



2 v^+25:*^'" 

Next we construct Eqs. ( 2.12D , ( |2.13| ), ( ^.14 ) as Hamilton equations by imposing the following Dirac 
brackets^ on independent variables at LF: 

(2.21) 
(2.22) 
(2.23) 
(2.24) 
(2.25) 
(2.26) 

while all other brackets vanish. The structure of these brackets evidently indicates that our indepen- 
dent canonical variables are not independent modes yet. Therefore as independent modes we propose 
to take the following linear but nonlocal combinations of fields: 

Q = (2.27) 

m 

Tfi 1 

= -r J * {IT - 2A + ^^Bi) A , (2.28) 

A_L — m 

Ci = Bi-di—^ -2\ + djBj) . (2.29) 

A_|_ — 

They satisfy the equations of motion 

(29+a_-A^ + m2)(^ = -^^±^^*j+-l(5_r , (2.30) 

(2a+a_-Ax + am2)Q = ^^^^-^ * j+ + 1 (a_r + 9./) , (2.31) 

m 2o- m ^ ' 

(25+5_ - A^ + m2)C, = f + d~ * j+ , (2.32) 



We have used relations ( p.lq ) and ( p.l6[ ) to remove dependent variables from the Hamiltonian 
^Though we have not used Dirac's procedi 
straints, therefore we call them Dirac brackets 



^Though we have not used Dirac's procedure of quantization, these brackets are evidently consistent with all con- 
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and their Dirac brackets have a diagonal form 



2d^_{C,{x),C,mDB 
2d-_{Q{x),Q{y)}^s 



-Sij6^{x-y) , 
-6\x-y) , 
6^{x-y). 



(2.33) 
(2.34) 
(2.35) 



The Hamiltonian density, when expressed in terms of the independent modes, looks like 



n 



eff 
jcov 



1 — a 



mQ + 



1 



25_ 



■J 



1 



1 



+ i;{diCi + — *j 



m? ^r, 1 r, m? n fjD + Q 
^(diQ) -—Q +i;{diV) + 



A 



m 



m" „<, 1 
2"" ' 2' 



d-r + dif) -fCi 



m 



■Q + 



m 



2d. 



■3 



(2.36) 



and, in spite of the presence of terms quadratic in currents j^, it describes the same physical system 
as the primary Lagrangian with constraints. This allows us to believe that the canonical structure for 
independent modes of vector fields, that we have just formulated, will survive in the complete theory 
with fermion fields. In order to incorporate our hitherto obtained results into the full interacting 
theory, we construct the effective Lagrangian density by means of another Legendre transformation 



pjmass 



d-dd+d - d+Qd-Q + d^ipd+^ - H 



eff 
jcov 

1 — a 



mQ + 



1 



■3 



1 



A 



m 



m 



Q + 



Ai + 



m 



-V — * 7 



2d 

1 

V + Q 



1 



:,{diCi + — *j 



d. 



2d- 



m 



d-j + dif 



(2.37) 



One can easily check that the brackets ( |2.35 - 2^35| ) and the equations of motion ( 2.30| - 2^3^ ) will follow 
directly from this effective Lagrangian within the canonical procedure without constraint^ 



2.2 Interactions with fermions 

Our next step, from the effective Lagrangian for the gauge sector to the full interacting theory with 
fermions, can be done easily by inserting fermion currents instead of and addition of the fermion 
kinetic terms 



= d_Cid+Ci - d+Qd^Q + d.ipd+^p + iV2iJjld+ip+ + V2ijj. 



id- — ed- 



(f + Q 



m 



V'- 



1 — a 



mQ 



2d- 



* r 



2 V 5_ 



+ 



A 



m 



m 



Q 



A_L +m^ 



m 



1 



2d. 



(2.38) 



^Strictly speaking there would be trivial primary constraints which usuaUy appear for covariant relativistic fields at 



LF. 
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where 



Ci + di 



^ + Q 



m 

Ci + a 



ni 



(2.39) 

(2.40) 
(2.41) 



As usually in the LF formulation, -!/;_ and fermions are non-dynamical and can be expressed in 



terms of dynamical fermions ip+ s-^id if). 



^ = J_ \ 



V2 ^^_ - ed. (^) 



(2.42) 
(2.43) 



and one obtains the nonlocal Hamiltonian density expressed solely in terms of dynamical independent 
modes 



'Htotal 



1 — a 



mQ 



1 



29_ 



1 

1 



* ,r 



2 2 2 2 V2 iO„ - ed- 



Ax -m2 Ax + m2 1 , 

Q + ^ + 

m m O- 



1 



2d- 



J+. 



(2.44) 



As we have expected, the non-vanishing equal (anti) commutators have the forms compatible with 
the former vector field brackets 



2d'^[^{x),^{y)] = -i5\x-y) , 
2&t[Q{x),Qm = i5\x-y) , 



{^:^(f),V+(y)} 



V2 



K+5^{x-y) 



(2.45) 
(2.46) 
(2.47) 

(2.48) 



and, just like in the ET approach to the covariant gauge condition, there are negative metric ex- 
citations here.^^ From these relations one can derive the dynamical equations for the interacting 
system^ 



(2a+5_ - Ax + m^)^ = 

Im 
^/2e 



1 



2m 



di 



* C + C * 



id_-e( ^) id--e( ^) 

\ m J \ m J 



Opposite signs in (2.46) and (2.47) allow the combination of fields if + Q to commute with itself on LF and this 
simplifies considerably the consistent definition of the integral operator [id- — ed- ("^^^)]~^- 

^^Here we do not discuss explicitly the ordering problems and the proper definition of singular products for noncom- 
muting operators. These very important aspects of the complete definition of Quantum Field Theory remain out of the 
scope of this work where we are ultimately interested in the Feynman rules for perturbative calculations. 
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A_L + 1 

m 2d- 



(2.49) 



{2d+d- - Ax + am^)Q 



V2e 
2m 

V2e 



9- 



2m 



di 



^d- - ^ (^) ^d- - e (^) 



* ^ + c * 



id- 



(^) 



(2a+a_ - Ax + m2)Q 



A_L - am^ 1 

2al 





m 




2m 








- * 


1 




72 





* J+ , 



(2.50) 



z9_ 



id- 



\ m J 



-idia' + M(3 + ea'd + a'di 



(f + Q 



m 



eV2 1 

+ —^7^ * 
2 5_ 



- ed- (^) 

^"-""^'Q + ^^v^ + a.a-(i-a)^*j^ 



(2.51) 



(2.52) 



V2 id_ - ed. (^) 



iV2 ,+ 1 



+ 2 "^^a 



A I — am^ 



i di a' + M/3 + ea^Ci + a'd, 
A_L + 



V^ + Q 



m 



m 



Q + 



m 



ip + dia-{l-a)—*J+ 



These equations show that one can introduce another quantum field A 

A = mQ + e-^ * J+ , 
2d- 

which satisfies the free field equation 

{2d+d- - A_L + am2)A = , 
but at the price of having extra non-vanishing commutators on LF 



il>\{x),K{y)\ = -e5^{x-y)^\{x) 



(2.53) 
(2.54) 

(2.55) 

(2.56) 
(2.57) 



This alternative, of having either the non-free field Q which commutes with fermions or the free field 
A which satisfies ( ^.56 ) and ( 2.571 ), is a specific feature of the LF formulation - contrary to the ET 
case where A is the only opportunity. 

The mere presence of q-number commutators would destroy a simple picture of perturbative calcu- 
lations based on Feynman diagrams and Wick contractions, therefore in the next subsection we will 
build the perturbation theory taking Q as an independent field and going to the interaction picture, 
where only c-number terms will appear as contractions for all fields. 
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2.3 Perturbation theory 

The perturbative calculations of the S-matrix elements are most easily performed in the interaction 
picture where all quantum field operators have free dynamics while interactions appear in the evolution 
of quantum states. Below we will work in this representation; however for clarity we will omit 
subscripts i for the field operators. Here the interaction representation is defined by taking the free 
Hamiltonian Hq as the limit of the total Hamiltonian ( |2.44 ) 

no = ]hnntotai = \{diCif + -^{diC,-djCif + "^0^ -]^{d^Qf -a'^Q^ 

+ 1(5,(^)2 + !^^2^_L^t_i_,^^^ (2.58) 

where 

eo = (-idia' + MI3) ^+ , (2.59) 
= {idiiP^a' + MiPl_(3) . (2.60) 

Next the interaction Hamiltonian, which is defined as the difference of these two Hamiltonians, can 
be written, for the later convenience, as a sum of two contributions 

T~(-int = 'Htotal — T~iint = T~^]nt + '^fnf (2.61) 

The first part describes the interaction of fermion and LF spatial components of the vector field 



+ -^4a^ (C, + djcl^)-^-^^ * a\C^ + 9,<A)V+ , (2.62) 

while the second one is connected with the current 

= -(2a,Ci + my.-amg + Ax0)^* J+-i^J+^* J+, (2.63) 



where we have introduced another notation 

, ^ + Q 



m 



(2.64) 



2.3.1 Field operators in the interaction representation 



From the free Hamiltonian ( |2.58| ) and the equal-x"*" (anti)commutation relations one derives the free 
field equations 

(2.65) 
(2.66) 
(2.67) 
(2.68) 
(2.69) 



(29+ 5_ - 


- A_L + m'^)Ci 


= 0, 


(29+ 9_ 


-A± + m'^)(p 


= 0, 


(2a+9_ - 


A_L + am^)Q 


= 0, 


{2d+d^ - 


Ax + M2)?/;+ 


= 0, 


{2d+d- - 


Ax + M'^)^l 


= 0. 
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Such free fields have their Fourier representations for all and for vector field mode^ we write 



aix) 

Q{x) 

Lp{x) 



(27r)2 Jo Ank. 



-ik-x 



k+-- 



oo ^2 



d^ki 



_oo (2vr)2 Jq Auk 



dk. 



-ik-x 



q{k) + e+'^-''q\k) 



k+-- 



d^k\ 



dk^ 



27r)2 Jq Airk. 



-ik-x 



p{k) + e+'^-'^p^k) 



where the creation and annihilation operators have non-vanishing commutators 



qik),qHk'] 
c,{k),c\{k') 



p{k),p\k') = {2^f 2k- &\k - k') 



(2^)3 2A:_ 5ij5^{k - k'). 



(2.70) 
(2.71) 
(2.72) 



(2.73) 
(2.74) 



Now it is an easy exercise to calculate the chronological (in x+) products of these free field operators 

(0|r+C,(x)C,-(y)|0) = 5ij^F{x-y,m') , (2.75) 
(0|r+<^(x)(^(y)|0) = A^(x-y,m2), (2.76) 
(0|r+Q(x)Q(y)|0) = -Ai.(x-y,am2) , (2.77) 



where the covariant Feynman massive propagator function /S.f{x,^'^) is defined in Appendix |B.l| . In 
the interaction Hamiltonian, the linear combinations of independent vector modes 



1 



* [2diCi + m(p — amQ + A±(f)] 



2a 

5-0, 

Ci{x) + di(j) , 



(2.78) 

(2.79) 
(2.80) 



are coupled with the fermion currents, therefore in the Dyson- Wick perturbation procedure we effec- 
tively encounter contractions given by the chronological products of B^^. Now after some algebra we 
find 



{0\TB^{x)B,iy)\0) 



-g^,Ap{x-y,m^)+d'py 
1 - a 1 



Af{x — y, m?) — Af{x — y, am? 



* S{x - y) , 



4 (a_)2 

where the last non-covariant contribution arises when one uses the equation 

{2d+d^ - Ax + y?)AF{x,p'^) = -i5^{x) 



(2.81) 



(2.82) 



in transforming the result into the form with second derivatives 5+ of the propagator functions. 
In later discussion it will be very convenient to reintroduce the dependent fermion fields 



1 1 



^/2 



id- 



(2.83) 
(2.84) 



^^The case of free fermion fields was given by Yan tlierefore we will omit them here and only quote all the needed 
results. 
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and consider the complete spinor fields = + ^- and = + ^1 . Taking the well known 
chronological product for independent fermions [|lO|] 



iV2A+dlAF{x -y,M'^ 



we derive the chronological products for dependent fermion fields 



+ 



rV'+(x)V'l(y) O) = UdfY + M)j°AF{x-y,M^) 



r+V-(x)vUy) 

(O |r+V-(x)VUy)| O) = ia^7+7°A^(x - y, M^) - l^l^^i^ * S{x - y) , 
and finally for the complete fermion fields 

(o|r+V'(x)V't(y)|o) = (idl^^' + M) jOAF{x-y,M^)-j+j^^*6ix-y) , 



or 



{0\Ti,{x)i;{y)\0) = {id^^j'' + m) A^(x - y, M^) - 7+-I- * 5(x - y) 



(2.85) 

(2.86) 
(2.87) 

(2.88) 
(2.89) 



Now we may re-express the interaction Hamiltonian in terms of B^, ip and ip^ . First we take Tij^^ 

(2.90) 



and, using the identity 
1 



id- — eB. 



1 

id- 



id- - 



-eB- 1 * ^ = ^ * ( eB^ 



we write it as 



int 



1 



id- - eB- 
1 



id- — eB- 



eB- ) * — * ^0 



V2^^' * id^ 



1 



* 1 + 



id- - eB- 
1 



id- — eB- 



eB- ]* — *^Q 



eB- * a'Biip^ 



(2.91) 



1 



+ -^^l^la'Bi*- 

V2 id- — eB- 



Then and can be expressed in terms of V'- and ^pl_ fields and finally, we arrive at the factorized 
form 

e , 1 



HL = i^e {j-B-+YB, 



1+ 



2 ' id-- eB- 

The second part of the interaction Hamiltonian is much simpler 

1 — Q ^ I 1 



7 B- +J^Br 



* 1p. 



'Hint = ^pe-f+B+ip 



(2.92) 



(2.93) 



where for simplicity, in the last term we have left the notation J+ for fermion current. 
These formulas show that here the noncovariant terms appear both in the propagators and the 
interaction Hamiltonians, contrary to the equal-time results where all the corresponding expressions 
are explicitly covariant. However, following the analysis by Yan |10[, one can hope that also here 



these non-covariant terms will cancel in pairs ( p.63| ) with ( 2.81 ) and ( |2.8^ ) with ( 2.92 ). 
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2.3.2 LF perturbative calculations of the S-matrix elements 



Here we will check the above conjuncture by studying the formal structure of perturbative calculations. 
The functional techniques [^] , |jl^ are quite useful for this purpose because they allow us to analyse 
contractions of gauge and fermion fields separately. First we check the contractions of vector field 
i?+ and treat all other components of vector fields and the fermion fields as classical objects for time 
being. The Wick theorem for transformation of chronological products into the normal products says 



T+exp 



where 



1 — a 



•+ 1 - 



exp 



exp 



(2.94) 



P++(x) = i (0 \TB+{x)B+{y)\ 0) + * S{x) 



ididl 



^f{x — y, m?) — Air(x — y, am? 



(2.95) 

and as usually, colons denote the normal product. Therefore we can simultaneously omit the instan- 
taneous current interaction in Hf^^ and take the covariant propagator for contractions of all vector 
field components 



-ik-{x) 



(27r)4 2k+k^ -k\-m? + ie 



-gtiu + (1 - a) 



kfj,ki/ 



2k^k^ 



/c^ — om^ + ie 



Thus effectively the complete interaction Hamiltonian is bilinear in fermion fields 

1 



H., 



int 



-0 * ej^Bfj 



2 lO- — eij- 



* ^ = -0 * V[-B^] * ip. 



,(2.96) 



(2.97) 



Next when all the vector fields are kept as c-numbers 6^, one can easily study the contractions of 
fermion fields 

(2.98) 



r+ exp -iV * V[b^] * V = exp [Tr In {1 - Sp * V)] : exp l-iip *V*{1-Sf* V)~^ * ip 
where now 



iSpix.y) 
iSpix) 



1 



(0 \Tip{x)ip{y)\Q) = iSpix - y) - * ^{x - y) 



One can check that the following factorization property holds 



(2.99) 
(2.100) 



1-Sf*V 



1 



Sf 



7' 



(1 - SFe^b^) 



7+ 



1 + 



1 + 



1 



7" 



2 id. 



eb. 



* e^'^bp 



2 id- - eb. 



* eYbp 



{1-Sf* ej^'b^) * V2 , 



(2.101) 



and then it is easy to notice that the non-covariant factor H2 disappears from the normal product 
part of Eq.(|]9|). It still formally remains in the closed loop contribution 



expTr In 1 + 



1 



2 id- - eb. 



expTr In 1 + 



1 



id- — eb. 



* eb- 
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However, the last expression can be shown to be independent of 5_ .[^ In this manner we have shown 
that the formal perturbative series based on the noncovariant interaction Hamiltonian ( 2.98 ) and the 



canonical noncovariant fermion propagator (2.99) will give the same result as the calculation based 
on the covariant interaction Hamiltonian ^ * e^^b^ * ip and the covariant fermion propagator Sp 

T+ exp -ii} * V[h^] * V = exp [Tr In {l-Sp* 67^6^)] : exp *{1-Sf* eYK)~^ * V'l : (2.102) 
provided all divergent expressions are properly regularized. 



3 LF Weyl gauge 

Another gauge condition which we have chosen for the massive vector fields is the LF Weyl gauge 
Bj^ = B~ = 0. This gauge condition is particularly suitable for the canonical procedure because it 
explicitly removes this component of vector field whose canonically conjugated momentum is zero. 
Contrary to the previous covariant gauge, the LF Weyl gauge can be strongly implemented in the 
Lagrangian density thus reducing the number of field variables. 



3.1 Vector fields sector 

First we take the model of massive vector fields coupled to the external currents and impose explicitly 
the LF-Weyl gauge condition. We see that the covariant mass term takes the form — ^m^(-Bi)^ 

^Tw%i = 0+B,{d.B, - d,B.) + i {d+B.f - i {d,B, - d,B,f - ^Bf + B.r + B,f (3.1) 

and all Euler-Lagrange equations are dynamical 

{2d+d^ - Ax + B, = di{d+B_ - diBi) + f , (3.2) 
d+{d+B^ - diBi) = j-. (3.3) 

Thus if we rewrite these equations in the first-order form 

{2d+d--A^)Bi = d,U + f, (3.4) 
d+U = j- , (3.5) 
d+B- = n + diB, , (3.6) 

and find the canonical Hamiltonian]^ 

,2 



n~ = U-d+B. + U%B, - C]^^%i = i(n)2 + ^{diB.f + mB. + ^Bf - Bd' - B^r, (3.7) 



^''One can proved perturbatively showing that the closed loop diagrams disappear, but in the present model, from 
6_ — d-(j}, one also gets 

[id- — e{d-(f>)]~^ {x,y) — expie(j>(x)(id-)~^ {x,y) exp—ieifiiy) 

and its functional determinant is evidently independent of 0. 
'^^The canonical momenta are simple here 

n- = d+B- - d,B, , 
= d-B,~d,B-. 
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then the Poisson brackets will follow immediately 



{2d.Bi{x+,x),Bj{x+,y)]p^ 
{S_(x+,x),n(x+,y)}p^ 



5ij5^{x-y) , 



5^{x-y). 



(3.8) 
(3.9) 



Just as in the previous case, our independent canonical variables {Bi, B^) are not independent modes. 
Rather we should take their linear combinations 



Ci 



Bi + di 



1 



A, 



* n 



A_L — m? 

which satisfy separated equations of motion]^ 
{2d+d^ - Ax + m^) d 

(9+n 



2djBj + Aj 



A_L — 



* n 



m? — Aj 
1 



* d-i , 



m 



A 



A_|_ — 
1 



* n 



Ai 



■A 



and still have diagonal Poisson brackets 



{2d^Ciix+,x),Cj{x+,y)}p^ 
{C-_(x+,x),n(x+,y)}p^ 



-6ij6^{x -y) , 
S^{x-y). 

At last we can express the Hamiltonian density ( |3.7D in terms of independent modes 

,2 ^2 



(3.10) 
(3.11) 

(3.12) 
(3.13) 

(3.14) 



(3.15) 
(3.16) 



njmass 

lit 



^total 



1 

2 

n- 



m"„ 1 

— n- — 

2 A| - 



* n - c_j" 



1 



A, 



dif + Aj 



1 



A 



* d-f 



1 



A, 



(3.17) 



and notice that no instantaneous interaction of currents occurs in this model. We see that the fermion 
field contribution is quite similar to that discussed in the previous case of covariant gauge and no 
modification of the perturbative vector field propagators should appear here. Therefore all we need 
to know here are the free vector field propagators. 

3.2 Free quantum fields 

We restrict our discussion to the free field case because in the interaction representation field operators 
have free dynamics. From the canonical analysis we take the form of canonical commutators at LF 



[2a_ C, (x) , Cj {y)] = -i5,,5^ {x - y) , 
[C^{x),U{y)] = i5\x-y) , 



and the free Hamiltonian density 



1 



(3.18) 
(3.19) 



(3.20) 



^The field C_ is a noncovariant multipole field which is characteristic for noncovariant gauge conditions. 
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They generate free field equations 



(28+ 



d+U 



, 
, 



1 



-m 



* n. 



(3.21) 
(3.22) 

(3.23) 



A_L — m? 

All fields, the covariant Cj and the noncovariant (11, C^) have their Fourier representations for all x 



Ci{x) 

n(f) 



-oo (27r)2 Jo Airk 



-ik-x 



■■Ci{k)+e+"'-''clik) 



d^ki 



dk- 
'2^ 



-ik-x 



p{k) + e+'^-^p\k) 



+ 



°° d^k\ 



(2vr)2 



- 

°^ dk^ 



n(x) 



27r 



ik-x 



c.{k)+e+'^-^cl{k) 



and the commutator relations for the creation and annihilation operators are 



c-{k),p\k') 
'ciik),c]ik'] 



cUk),p{k')\ = i{2TTf6^{k - k') 
{2nf 2k^ 6ij5^{k-k'). 



(3.24) 
(3.25) 

(3.26) 

(3.27) 
(3.28) 



Now a quite straightforward calculation gives the chronological products for independent modes 

(3.29) 
(3.30) 

(3.31) 



< 0|r Ci{x)Cj{y)\0 > = dijApix - y) , 
<0|TC_(x)n(y)|0> = El{xL-yL)5\x^-y^) 



< 0\T C_(x)C_(y)|0 > 



m'^Ep{xL - yh) 



1 



Ai 



* 5 (x_L - y±), 



where all other propagators vanish^'^. We notice that the causal propagation of noncovariant fields (H, 
C- ) takes place only in the LF longitudinal directions xl = (a;+ , x^ ) and in the transverse directions 
x_L = {x2.,x^) these propagators are either local ( |3.30| ) or are given by the modified inverse Laplace 
operator (|3.3l| ). This means that their causal properties i.e. the ML-prescription which appears in 
them, do not depend on the space-time dimensionality, and thus have no direct infrared singularities 

Finally, in order to find the propagators for primary vector fields, we use the relations 

1 



Ci-di 



A 



* n 



- 

1 



Ai 



2a,- C,- - A^ — 
^ ^ A 



1 



*n 



and then after some algebra we obtain the expressions 



< 0|r Bi{x)Bj{x)\0 > 



SijAp^x -y,m'^) , 



(3.32) 
(3.33) 

(3.34) 



B.l 



^^Definitions and properties of tiie noncovariant propagator functions Ep (x) are given in Appendix 
^® Evidently tlie infrared singularity would appear in the integral operator (Ax — m^)~^ when the limit ^ is 
taken in two transverse dimensions. 
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<0\T B^{x)Bi{y)\0> 



< 0\T S_(x)B_(y)|0 > 



df- 



Ai 



df 
29^ A 



x^—y 







2dlAF{x -y) + El{xL - yL)5\xi_ - y^] 

d^Apitx-y) , 
1 



(3.35) 



(Aj 



m 



2\2 



m Ep{xL - yi) 



1 



A_L — rri 



2dlAF{x -y) + Ej,{xL - yL)S\x^_ - y^] 
* S'^{x± - y±) 



29^ 



x^-y 



d^AFi^,x-y)+m^ 



x^-y 



dS, / d-qApi'n.x-y) , 



(3.36) 



with the successively increasing number of independent modes contributions^^ These all components 
have a concise form in their Fourier representation 



< 0|r Bf,{x)B,{y)\0 > 



d^k 



Jk-(x-y) 



+m 



(27r)4 A;2 - m2 + ie 
+ ie'sgn(A:_)]2 



k^ + ie'sgn(A;_^ 



(3.37) 



where we have introduced the LF Weyl gauge vector A'^ = (A^_ = 1,A^+ = ^■,N±_ = 0) which 
chooses the LF-Weyl gauge condition N^B^ = B^ = 0. In this way we have obtained the causal 
ML-prescription for the spurious poles in the vector field propagator via the canonical quantization 
procedure on a single LF surface of quantization. We see that the noncovariant modes 11 and C_ are 
inevitable for this encouraging result. Thus we can speculate that the canonical procedure at a single 
LF cannot lead to the ML-prescription for the LC-gauge, where all noncovariant nonphysical modes 
are excluded. 



In the second and the third case we have used the property (B.9) of the propagator functions Af{x) i Ep{x). 
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Part III 

Light Front QED 



4 Class of LF Weyl gauges 

Now we begin the discussion of the Abehan gauge field models where the gauge transformation is a 
local symmetry of the theory. Therefore now it is not a problem of having smoother behaviour of 
canonical LF commutators but the fundamental property of a gauge field prescription which makes us 
choose some gauge fixing condition. The first choice that we discuss here is the class of LF Weyl gauge 
conditions which is a generalization of the exact LF Weyl gauge = 0. It is usually introduced into 
the Lagrangian density by the so-called gauge fixing term 



1 



— {N^^A.f, (4.1) 

where in our case the LF Weyl gauge vector is A^^ = (A^'*" = 1,A^~ = 0,A^-'" = 0). This term is 
very convenient in the path-integral approach to the quantization of gauge field theories |32] because 
it constitutes a non-singular quadratic part of gauge field Lagrangian thus allowing for immediate 
inversion. However the path-integral procedure gives no prescription for spurious poles in the gauge 
field propagator which should have the form 



D^uix) = i{0\TA^{x)A,{0)\0) 

d'^^ -ik-x I 1 



(27r)4 1 2k+k- -k^ + ie 



[k+ + ie'sgn{k-)] 



kit 

a- 



[k+ + ie'sgn(/c_)]^ j 

(4.2) 

with the causal ML-prescription + ie'sgn(A;+)]~^'^. We see that while the above propagator is 
finite for a ^ where it describes the exact LF Weyl gauge, the expression ( [4.1D is ill-defined for 
a = 0. Therefore we choose another form of the gauge fixing condition via the Lagrange multiplier 
field A 

AA+ - |a2, (4.3) 

which is equivalent to the previous choice as long as a 7^ 0. Evidently it is also regular for a = 
when it describes the exact LF-Weyl gauge A^ = 0. In the present case A is not a dynamical field, 
contrary to the case of covariant gauges. Below we will consecutively discuss the sector of gauge fields 
in two distinct cases of 1+1 and D+1 dimensions. 



4.1 Gauge field sector in 1+1 dimensions 

We want to discuss first the case of 1+1 dimensions where the infrared singularities in the transverse 
momenta for A:^^ = are excluded from the considerations. Thus we can focus our attention on 
singularities in the longitudinal momentum Just as we did for previous models, also here we start 
from the sector of gauge fields coupled linearly to external currents 

^Iweyi = ^(9+A_-a_A+)2 + A_r + ^+i+ + A^+-fA', (4.4) 

and the Euler-Lagrange equations are 

d+ {d+A^ - d-A+) = j- , (4.5) 
d- {d+A. -d^A+) = -j+ - A , (4.6) 
A+ = qA. (4.7) 
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Without going into details, we give the canonical Hamiltonian structure which follows from the above 
equations after removing non-dynamical field variables (yl+ , A) . The canonical Hamiltonian density 



n 



can, 1+1 



hi- (i - adl) n- - A^r + aj+9_n- + 



aWeyl c^' 

and non-vanishing commutator at LF 

[n(x+,x"),yl_(x+,y")] = -i 5{x'' - y'). 
generate effective dynamical equations 

= J" 



[l-ad'i)ir -ad-j^. 



(4.8) 



(4.9) 



(4.10) 
(4.11) 



We stress that this quantum theory describes a larger system than the physical excitations. This is 
due to the lack of the Gauss law as an equation of motion 



G = d-{d+A^ - d-A+) + j+ = d-U- + f 



0. 



(4.12) 



Evidently one cannot strongly impose the condition G = because this would be in a conflict with 
the commutator ( ^.9[ ) ; thus one has to implement it weakly as a condition on states [ ^ ] 

{phys'\G{x)\phys) =0. (4.13) 

indicates that external currents are coupled with the 



The interaction part of the Hamiltonian 
linear combinations of fields 



A_ 
A, 



A_ 



-ad-U- 



(4.14) 
(4.15) 



and also here there is an instantaneous interaction term ^j^j^ which, during the Wick contraction 
procedure, would modify the perturbative propagator from the canonical chronological product to 
the following expression: 



D^,{x -y)= i(0|r+ i^(x)i,(y)|0) - ag^^g^j'^ix - y). 



(4.16) 



First we calculate the chronological product for independent modes in the free field case when they 
can be represented by 



ir{x) 
A_{x) 

with the Fourier representation 



vr X 



x+{l - ad'i)-K{x~) + a_(x") 



a-[x 



vr X 



°° dk_ 







27r 



-ik—x 



-ik—x 



a{k-) + e 
p{k-) + e 



a^{k^ 



+ik—x t 



p*{k^ 



and the commutation relations for the creation and annihilation operators^ 

'a-{k_),p\k'_)\ = \a^_{k^),p{k'_)] = i 2tt 5{k_ - k'_) 



(4.17) 
(4.18) 



(4.19) 
(4.20) 

(4.21) 



^"The weak Gauss law condition ( [4.13| ) shows that for free fields, all physical states are created by the p(fc)^ operators 
and hence they have zero norm. This means that no physical photons are present in 1+1 dimensions. 
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Now the free propagators are 

(0|T+ n-(x)n-(y)|0) = 0, (4.22) 
{0\T+ A^{x)U~{y)\0) = E^x-y), (4.23) 
{0\T+ A^{x)A^{y)\0) = {I - adl)El{x - y), (4.24) 

where noncovariant functions Ep{x) and Ep{x) are defined in Appendix |B.1| . The perturbative gauge 
field propagator has the causal form 

and simultaneously the current interaction term is omitted in the interaction Hamiltonian. Then one 
could incorporate the fermion interaction, but we will not discuss it here because this would be a mere 
repetition of the relevant considerations from the subsection ( |2.3.2| ) with a rather evident neglect of 
transverse coordinates and components. 

We conclude that in 1+1 dimensions, the class of the LF Weyl gauges effectively leads to the causal 
form of the perturbative Feynman rules though at the canonical level it contains non-covariant terms 
which ultimately cancel in pairs. 



\k+ + ie'sgn(A;_)l 



a 



\kA. + ie'sgn(A;_)P 



(4.25) 



4.2 Gauge field sector in D+1 dimensions 

As a physically relevant model we take the higher-dimensional case where excitations of physical 
photons are possible. However if we would take 3+1 dimensions then we would encounter infrared 
singularities connected with the inverse Laplace operator in 2 transverse directions. Therefore we have 
decided to use here the dimensional regularization of this singularity by working with d = D — 1 > 2 
transverse coordinates x±_ = (x2, a^s, . . . , xd). Thus we would like to start with the Lagrangian density 

Cc^weyi = {d+Ai-diA+){d.A,-diA_) + ]^{d+A^-d-A+f --^{diAj-djAif 

+ A+A-^K^ + A+j+ + A.r +Ad\ (4.26) 

which generates the Euler-Lagrange equations 

d+ {d+A^ - d^A+ - djAj) = -A;[A+ + j" , (4.27) 
-d^{d+A_ - d^A+ + djAj) = -A;[A_ + i+ + A , (4.28) 

{2d+d- - Al) Ai = di {d+A^ + d-A+ - djAj) + f , (4.29) 

A+ = aA , (4.30) 

where A^ = {dj)"^ denotes the Laplace operator in d > 2 dimensions. As the independent modes we 
take the modified canonical momentum 

n = d+A_ - diAi - d-A+ (4.31) 

and the gauge fields 



Ci = Ai-di^*{U + 2djAj) , (4.32) 
C- = A^-d-^*{U + 2djAj) (4.33) 
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thus the dynamical equations of motion are 



1 



1 



The canonical Hamiltonian density can be expressed in terms of the independent modes 



(4.34) 

(4.35) 
(4.36) 



n 



eff 
aWeyl 



a 



n-i- * 



- Ci 



1 



f-2di^*{d-r + d,f) 



dif + 

c-j- , 



and the non-vanishing commutators at LF are 

[n(x+,x),C_(x+,y)] 
[2d.Ci{x+,x),Cj{x+,y)] 



^- ■ 



i6i,5''^\x-y). 



(4.37) 



(4.38) 
(4.39) 



Just like in the low-dimensional case, we have the instantaneous current interaction which, during the 
Wick contractions, modifies the perturbative gauge field propagator from the chronological product 
form 



D^,,{x -y)= i(0|r+ A^{x)A,{ym - ag-^g-,5''+^{x - y) 
where we have introduced the notation 

Ai = Ci-d~*{Ii + 2djCj) , 

1 



(4.40) 



C--d-^*{Ii + 2djC,) , 



(4.41) 

(4.42) 
(4.43) 



A- 

A+ = -aAfC-. 

The free propagators for independent modes are calculated directly from the Fourier representations 



C_(x) 

a{x) 

n(x) 

7r(x) 



c_(x) = 



oo 



oo {27ry Jo 2-jr 
°° dk- 



ik-x 



a{k) + e+'''-^a\k) 



J {2-kY Jo 2tt 2k- 
7r(f) -ax+(Al)2c_(f) 



ik-x 



^^c,{k) + e+"'-^4{k) 



k+-- 



d'^ki 



dk. 



-ik-x 



p{k)+e+'^-^p'^{k) 



/-oo {2-kY Jo 27T 

where the non-vanishing commutators for the creation and annihilation operators arep^ 



a{k),p^{k" 



a^k),p{k')\ =i{2^f+^5'^+^{k-k') , 



Ci{k),c]{k') = {2TTf+^ 2k., di^S^+'ik - k') 



(4.44) 

(4.45) 

(4.46) 
(4.47) 

(4.48) 
(4.49) 



^^We recognize that the present operators a' and p' are trivial generahzations of the respective operators in 1+1 
dimensions. Free field physical states, selected by the weak Gauss law, can be created by cj and excitations, where 
the former would be the positive norm photon states while the latter would be the accompanying zero norm states. 
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where 6'^'^^{k) = 6'^{k±)5{k-). Determination of free propagators for independent modes is rather 
simple 



(o|rc_(x)n(y)|o) 
(o|r n(x)n(y)|o) 

{0\T Q{x)Cj{ym 



El{xL-yL)5'^{x±-y±) , 
a{AifEUxL-yL)S''{x± 



y±) 



S,jDf^\x-y) 



(4.50) 
(4.51) 
(4.52) 



where the covariant Feynman propagator function D'p~'^{x) is defined in Appendix B.l, while the 
calculation of the perturbative gauge field propagator is quite tedious; here is the final result presented 
in concise Fourier representation 



Df,u{x -V) = j 



d'^+^k 



(2vr) 



d+2 



1 



{k^N^ + k^N^) 



/c+ + i e' sgn(A;_ 



a 



k^ky 



[/c+ + i e' sgn(A;_)]2 



(4.53) 



which evidently has a regular limit d ^ 2, where it gives the expected result ( [4.2| 



5 General axial gauge 

In this section we would like to discuss the general axial gauge condition imposed on the gauge field 
potential n^^A^ = 0, where the axial gauge vector has the form n"*" = l,n^ = — a, n"*- = 0. This 
general choice will allow us to analyse and compare within the LF canonical formalism different 
gauge conditions: 

• the LF-Weyl - for q = ; 

• the temporal Minkowski - for a = —1 ; 

• the spatial Minkowski - for a = 1. 

Also the verification whether the limit a — > ±00, can be considered as a possible limiting procedure 
leading to the LC gauge = 0. The present form of the gauge condition n^A^ = A^ — aA^ = can 
be implemented either explicitly or via the Lagrange multiplier field; below we take the first possibil- 
ity. Also here we discuss two cases in 1+1 dimensions and in 3+1 dimensions because they have quite 
different physical interpretations. Again we deal explicitly only with the gauge field sector because 
the discussion of the fermions and interactions with them would follow along the lines described in 
Section |^. 



5.1 Gauge fields in 1+1 dimensions 

Now the Lagrangian density is greatly simplified 

4-- = \id^A^-ad^A^f-+A^{j-+aj+), (5.1) 
and generates only one Euler-Lagrange equation 

{d+-ad.fA_ = j-+aj+ , (5.2) 
which is equivalent to the system of the first order equations 



d+U = ad-U + j~ + aj-^ , (5.3) 
d+A^ = U + ad-A.. (5.4) 
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When these equations are assumed as the Hamilton equations with the Hamiltonian density 



^con — n (9+A_ — Ctemp 

then the non-vanishing Poisson bracket is canonical and the quantum commutator is 

[U{x-),A^{y-)] =-i6{x--y-). 



^{U~f - aA^d^U- + -a{diA.)^ - A^{r + aj+) , (5.5) 



(5.6) 



Inspecting the Hamiltonian we find no instantaneous interactions of currents, thus the perturbative 
and free propagators will coincide here. The quantum free fields can be given as 



U{x) 
A^(x) 



7r{ax^ + X ) , 

a„(ax^ + x~) + x~^Tr{ax'^ + x~) , 



(5.7) 
(5.8) 



where the fields 7r(x) and a-{x) were defined by ( [4.19| ) and ( 4.20 ) in Subsection 41. Therefore the 
whole discussion given there applies also here and in order to avoid unnecessary repetitions, we just 
take the results 



(o|r^_(x)n(y)|o) 

(0|T^_(x).4_(y)|0) 



Ep[x^ — y , a{x~^ — y^) + x 
El[x+ 



y ,a^x 

Next we easily find the form of the perturbative propagators 



y ) + X 



y 



(5.9) 
(5.10) 



D__(x) = (0|r ^_(x).4_(0)|0) 
D+^{x) = a(0|T yl_(x)^_(0)|0) = i 
D++{x) = a^{Q\T A^{x)A^{{))\Q) = i 



ik-x 



-oo {2TTf 

°° (fk 



\k-\. — ak- + ie sgn {k-)Y ' 



-ik-x 



a 



» (27r)2 
°° (fk 



[kj^ — ak- + ie sgn (A;_)]^ 



ik-x 



a 



-oo (27r) 



[/c_l_ — ak- + ie sgn (/c-)]^ 



(5.11) 
(5.12) 
(5.13) 



which can be written concisely as the Fourier integral 



oo 



fk 



ik-x 



-n 



oo (27r)2 2k+k- + ie 
2 kfj^ky 



ak- + ie sgn (/c_ 



[k^ — ak- + ie sgn (k-)]'^ 



(5.14) 



We see that this propagator exists for all finite values of a with the ML-prescription for spurious 
poles. Thus even for the spatial axial gauge we can produce a gauge field propagator with the causal 
spurious poles when quantizing canonically at LF, contrary to the FT formalism where such case is 
not possible [^]. When one takes the limit a — > ±00, the causal nature of spurious poles in the above 
propagator is lost and one again has the LC-gauge case with the CPV poles. 
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5.2 Gauge fields in 3+1 dimensions 

Now we would like to analyse the general axial gauge condition in the more physically relevant case of 
3+1 dimensions. Here again the gauge condition is implemented explicitly in the Lagrangian density 

1 



Haipha = {d^A,-adiAJ){d.A,-diA^)^)^(d^A^-ad-AJj' --^{diAj-d.iAi, 



(5.15) 



where the external currents describe interactions with the charged matter. The canonical analysis 
starts with the Euler-Lagrange equations 



A. 



= {d+ + ad-) d^Ai + j + aj+ , 
i2d+d--A^)Ai = ^^[{^+ + a^-)A.-^JAj]+f 

which are equivalent to the Hamilton equations 

(9+n = ad- (n + 2diAi) - 2aA±A^ + j~ + qj+ 
(2a+a_ - Ax) Ai = diiU + 2ad.A_) + / , 
d+A- = U + d^Ai + ad-A- . 



(5.16) 
(5.17) 



(5.18) 
(5.19) 
(5.20) 



The canonical Hamiltonian density is 

Hcan = ^(n)2 + ^{diAjf + Ud^Ai - aA-d-{U + 2diAi) - a{diA. 



A.{j- + aj+) - Aif (5.21) 



with the Dirac brackets 



{n(f),A_(y)}^^ 
{2d.A{x),A,iy)}^^ 



-6ij6'^{x -y) , 



(5.22) 
(5.23) 



while other brackets vanish. When trying to separate these independent variables into the independent 
modes we encounter the problem of inverting the differential operator at LF A~ = 2ad'^ — A± which, 
only for a < 0, is an elliptic operator with the regular Green function (defined in Appendix B.3 ). 
Thus we have to choose only negative values of a which is equivalent to the selection of temporal 
axial gauges.]^ Now we define independent modes 



A = H - 2aa^ — * [H - 2ad-A- + 2djAj] 

Aq. 

C_ = A- + d-—*[U-2ad-A_ + 2djAj] , 

Aq 

Ci = Ai + di^*[U-2ad^A-+2djAj], 



which satisfy dynamical equations of motion 

id+-ad-)C- = ^* [5fe/ + a_(r +aj+)' 



{d+-ad-)A 
{2d+d. -Ai_)Ci 



2qA^ C-+j + aj^ 



2ad- * 

Aq 



5u' + 5-(j- + «i+) 



f + 2di—* 



9fc/ + a_(r + aj- 



(5.24) 
(5.25) 
(5.26) 

(5.27) 
(5.28) 
(5.29) 



^^The cases of spatial gauges [a > 0) will not be discussed here and are left for future investigations, while the null 
gauge (q = 0) is singular at 3+1 dimensions and needs some infrared regularizations. 



33 



and have non-vanishing commutators at LF 



2d-_[Ci{x),C,m 
[m,C-{y)] 



-i5ij5^{x-y) , 
-i6^{x — y). 



(5.30) 
(5.31) 



The Hamiltonian density ( 5.21] ), when expressed in term of these modes 



[A - 2ad-C- + 2djCj] 



(j 



C-+d-—*[A- 2ad-C- + 2djCj , ] 



(5.32) 



contains no direct interactions of currents, therefore the perturbative propagators are the chronological 
product of free fields. For free fields we write 



A{x) 
Qix) 



C-{ax~^ + X ,x±), 

2ax~^ A~c- (ax"*" + x~ ,x±) + X{ax~^ + x~ ,x±) 

Ci(x), 



where 



A(x) 



CiiX 



oo 



oo 



oo (27r)2 Jo 27r 
°° d^ki dk_ 



-oo (27r)2 Jo 27r 
°° d^ki dk. 



/-oo (27r)2 Jo 2tt 2A;_ 
with the commutators for creation and annihilation operators 



e-*^-^a(A?) + e+^^-^at(fc) 

e-*^-^p(A?) + e+^^-V(A?) 
-'''■''ciik) + clik) 



a{k),p\k') 
'c,{k),cUk') 



a\k),p{k')\ =i{2TTf6^{k-k') , 
(27r)3 2A;_ Sij5^{k-k') , 



while other commutators vanish. Now we easily find the relevant chronological products 

(0|rC_(x)A(y)|0) = EiF{^L-yL)S\x±-y±) , 
{0\T A(x)A(y)|0) = -2aA^ElpixL - yL)S\x^ - y^) , 



{0\T a{x)Cj{ym 



6ijDp{x - y), 



(5.33) 
(5.34) 
(5.35) 



(5.36) 
(5.37) 
(5.38) 



(5.39) 
(5.40) 



(5.41) 
(5.42) 
(5.43) 



where e]^p{xl) = Ey^{x'^,x + ax^). The perturbative propagators contain the following linear 
combinations of independent modes: 



n 

A. 
A, 



A - 2ad'i * [A - 2a5„C„ + 2djCj\ , 
C_ + — * [A - 2a5_C_ + 2djCj] , 
Ci + di^*\A- 2ad-C- + 29. a 



■A 



(5.44) 
(5.45) 
(5.46) 
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and we find after some algebra 

(0|r A_(x)^_(y)|0) 

{0\T A{x)A,{y)\0) 

{0\T Ai{x)Aj{y)\0) 

Next, using the relation A^ = a A. 
gauge field propagator 

-n2 



i2dld^^ (Df^eIp) {x-y) , 
i2dUad-_ + d+) {Df * EIp) (x-y), 
S^jDpix -y)+ i2adfd^ (Df * E^j, 



y)- 



(5.47) 
(5.48) 
(5.49) 



we find the Fourier representation for all components of the 



-ik-{x-y) I 



2k+k^ - k"^, +ie 



-9fMU + 



fc+ — ak- + i e' sgn(A;_ 



— ak- + i e' sgn(A;_)]2 



(5.50) 



Thus we have found the perturbative gauge field propagator with the ML-prescription for spurious 
poles for all temporal gauges. Now the LF Weyl gauge can be understood as the limit a ^ taken 
in ( 5.5C| ). One can compare all LF Weyl propagators obtained via different routines: the massive 



electrodynamics ( 3.37| ) when 



m 



0, the class of the LF Weyl gauges ( [4. 531 ) when a — > 0, and the 



above general axial gauge ( 5.5C| ) when a — > leading always to the same result 



{Q\T A^{x)A,{y)\{))LFWeyl = i 



d^k 



-ik-{x-y) 



(27r)4 2k+k- -k\ + ie 



-9iiv + 



{k^N^ + k^N^) 
k+ + i e' sgn(A;_) 



(5.51) 



The other passage limit to null gauge, when a 
CPV for the LC-gauge propagator. 



iboo changes again the ML-prescription into the 



6 Flow covariant gauge 

The analysis presented in Section ^ can be considered as an infrared regularized model which in the 
limit leads to the QED for the class of covariant Lorentz gauges c?^^'^ = aA. Therefore 

here, instead of repeating the previous analysis for explicitly massless vector gauge fields, we decided 
to study another class of gauge conditions: the flow covariant gauge d+A^ + ad-A+ — ad±A± = 
which for a = 1 produces the results for the Lorentz gauge while for a = it describes the LC-gauge 
A- = 00 

6.1 Model in 1-|-1 dimensions 

First we would like to discuss the electrodynamics with charged fermion fields in 1-1-1 dimensions 
which is classically described by the Lagrangian density 

J^'ftoT'' = \ (5+^- - d-A+f + ^ {ij^^d, - e^A^ - M) + A (9+A_ + a9„ A+) , (6.1) 

where the Lagrange multiplier field A implements the flow covariant gauge. Here we explicitly see 
that the flow gauge is attainable, by some suitable gauge transformation, for all values of the gauge 
parameter a ^ 1. Thus we expect that we may encounter expressions singular at a = — 1 during the 
canonical procedure and in final results. 



Strictly speaking it leads to the modified LC-gauge d+A- = which should have the same perturbative Feynman 
rules as the LC-gauge. 
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6.1.1 Gauge field sector 

According to our previous routine, we start with the sector of gauge field coupled with the arbitrary 
external sources 

C>^L = lid+A_-d^A+f + A^j- + A+j+ + A{d+A_ + ad^A+), (6.2) 

and here we have the Euler-Lagrange equations for the gauge potential and the Lagrange multiplier 
fields 

d+ {d+A^ - d-A+ + A) = j- , (6.3) 

d- {d-A+ - d+A_ + aA) = j+ , (6.4) 

d+A_ = -ad-A+. (6.5) 

However two fields can be parameterized 

by means of the canonical conjugate momentum field n~ and effectively there are only two dynamical 
equations of motion: 

a+n- = r , (6.8) 

with the canonical Poisson-Dirac bracket at LF 

{A^{x+,x-),U-'{x+,y-)}^s = S{x- -y-) , (6.10) 
while the canonical Hamiltonian density is 

^^r- = n"^-^+-4/^ = ^(ri^("n--^*j+)'-^-r. (6.11) 

Now we can give the effective Lagrangian density for the gauge field sector 

= n-5+^--2(ii^("n--^*J+)' + ^-r , (6.12) 
and then easily incorporate fermions back into the dynamical system. 

6.1.2 Interaction with fermion fields 

When the gauge fields are described by effective independent variables, the system containing complete 
dynamics of fermions and their interaction with gauge fields is given by the Lagrangian density 

^n^^"" = ^-d+A_ + V2i^pi^+^+ - j (an- + 6^2^- * (V^U+))' 

+ V2iipld^ip^ - eV2V'I^-A_ - M (V'l7°V'+ + V'+t"^-) • (6-13) 
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As it usually happens in the LF formalism, the fermion fields i tjj^_ are non-dynamical and as 
dependent variables can be expressed by other fields 



1 1 



y/2 id- - eA. 



(6.14) 
(6.15) 



In this way, we have arrived at the Hamiltonian density which depends solely on the independent 
dynamical fields (which are already independent modes) 



n 



1+lQED 
flow 



2( 



1 / 1 \ ^ M"^ 1 



(6.16) 



Now the canonical quantization is immediate, one takes the properly ordered expression (lel^ ) as the 
quantum Hamiltonian density and non- vanishing (anti) commutation relations 



[ll-{x+,x-),A^{x+,y-)\ = i6{x--y-), 
|V+(2;+,x"),V'+(x+,y")| = ^5{x'-y'). 



(6.17) 
(6.18) 



Here the equations for quantum operators will have the same functional form as their classical 
counterpart^^ 



d+A. 
-iV2d+i; 



. 1 1 

^"71^ * id- - eA_ id- - eA. 



a 



(l + a)2 
M2 



air + eV2^ * (^+V'+) 



1 



y/2 id- - eA. 



1 



t 



1 



^2 ^+ * 'id- - eA- (1 + a)2 ^+ d 



(l + a)2^+9_ 



aU- + eV2— * ('0+V'+) 
aH" + eV2^ * ('0+V'+) 



(6.19) 
(6.20) 
(6.21) 
(6.22) 



6.1.3 Perturbation theory 

In the interaction representation, the field operators have free dynamics generated by the free Hamil- 
tonian which in the present model is 



Hq 



l+lflow 



r2 

2(1 + a)^ ' ■ 

The remaining part of the Hamiltonian is taken as the interaction Hamiltonian 



(6.23) 



n 



I+IQED 
int 



n 



I+IQED _ y^l+lflow 



flow 







-^n--*(^U4 



(l + a)2 



M2 



V2 



1 



id- id- — eA- 



* V'4 



(6.24) 



^*The expression ?/)^(a:)i/)+(a;) is singular for quantum fields operators and needs some regularization. If 
one takes the splitting-point method which is compatible with the local gauge symmetry limr^-^o il>\.{x^ , x~ — 

-ie P_ '^"^ d^A_(x+ , _ , - i 

ri)e -"^ -1 +»?), then the equations for 11 , tf}+ i tp^ are modified. However, this would lead us 

beyond the scope of this paper and no such regularization will be discussed hereafter. 
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and we see that there is a direct instantaneous interaction of currents which will modify the pertur- 
bative propagators from the primary form of Wick's contractions of effective gauge field potentials 



V^,{x - y) = i (0 \TA^[x)A,{y)\ 0) + ^fip^^ * - v) 



where 



A. 
A, 



A. 



a 



1 



* n" 



We need the free propagator for independent modes and they can be obtained quite easily 

(o|r+ n-(x)^_(y)|0) = Elix-y) , 



a 



(0|T+ A_(x)^_(y)|0) = ———E%{x-y), 

(1 + ay 







iV2A+d'-lAl{x -y,M^ 



(6.25) 

(6.26) 
(6.27) 

(6.28) 
(6.29) 

(6.30) 



where A|,(x, M^) is given in Appendix [B.l| . Then one finds the perturbative propagator for the gauge 
fields 

rEUx) , (6.31) 



D__(x) 
D+-{x) 
D,,(x) 



a 



(l + a)2- 

Q 



1 



(1 + af d. 
1 1 



E^x) 
* 6{x — y) 



(6.32) 
(6.33) 



which has non-causal structure for (+) components. This is connected with the infrared singularities 
of covariant massless fields in the 1+1 dimensional LF formulation |^^. Specially the LC-gauge limit 
a — > produces the gauge field propagator with the CPV prescription for spurious pole and no causal 
pole at all. 



6.2 Higher-dimensional model 

The failure of the 1+1 dimensional approach towards perturbative propagators with causal poles 
indicates the danger of the infrared singularities also in the 3+1 dimensions. Therefore here we 
introduce the dimensional regularization by taking d = D — 1 > 2 transverse coordinates x± = 
{x2, . . . , Xd-i-i}- The gauge field sector described by the Lagrangian density 

jr^D+l gauge ^ ^d+A-d,A^){d^A,-d,A^) + ^{d+A^-d^A+f-^{d,A,-d,A,f 

+ A.r +A+j+ + A,f + A{d+A^ + ad^A+-adiA,), (6.34) 

still contains constraints; therefore, omitting all details which will be presented elsewhere, we can use 
the effective description 

^ d-ad+Q-^{diCjf-{l + a)d+Xd^<P + ad,Xd,<P 

- cp[d-j- + dif)+fCi 
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where the independent modes (Cj, (p, A) are defined by primary fields 



A 

C,; 



1 

A- 



(1 + a){diAi - d-A+) -aA- — *j- 



1 1 
1 + adI 
Ai - did). 



(6.36) 

(6.37) 
(6.38) 



Next we add directly the fermion contributions and as the starting point for the quantization of the 
complete model, we take the effective Lagrangian density 



rD+l 
'-'flow 



d-dd+d - (1 + a)d+Xd-(t> + iV2^pld+ip+ + V2ipl_ {id- - e5_0) ^p. 



a 



1 + aJ d. 



— * j+ - eV- - i^U , 



where 



^ = + M/3) V+ + e (Ci + 9i0) 

= (^idii;ia' + Mi;lf3)+e2pla'{Ci + di4 



After removing the dependent fermion fields ijj- and ijj. 



t 



1 



^/2 id- — ed- 



we have the Hamiltonian density 



njD+l 

'^flow 



^' V2 id--ed- 



1 



(6.39) 

(6.40) 
(6.41) 
(6.42) 

(6.43) 
(6.44) 



1 1 



(6.45) 



which depends only on the independent modes, and the non-vanishing (anti)commutators at LF are 

(6.46) 
(6.47) 

(6.48) 



2dl[C\{x),C,m 
{l + a)d-^[<P{x),m] 



-i6ij6^{x-y) , 
i6^{x-y) , 



V2 



A+5%x-y) 



Therefore we have the field equations 



(6.49) 
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[{l + a)d+d- -aA±](j) 
[(1 + a)d+d_ - aA±] X 



9 X a 



1 



V2 



di 



V2 
Aj 
1 

7i 



5_ 



j+ 



1 



(6.50) 



i5_ — e(?_ 
1 



1 



1 



a 1 
1 + adZ 



id- — ed-4> id- — ed- 
* J+ , 



id- — e9_ 



(6.51) 



-i^io^ + M/3 + e (Q + 9i(/)) a' 



1 



1 + 



* I aX + 



5_ 



1 1 

Q 



i5- 



39- 



1 + a 



(6.52) 



V2 
eV2 



idi a' + Mp + e{Ci + di(p) a' 



1 



id- — ed- 



1 + a'^^'5_ 



* aA + 



1 1 

1 + aaI 



a 



1 + a 



(6.53) 



which describe the quantum theory in the Heisenberg picture. We also notice that for the complete 
interacting theory one can define a free field 

1 1 



A = A 



1 + ad. 

which satisfies noncovariant dynamical equation 

[(1 + a)d+d- - aA±] A = 

but has two extra non-vanishing q-number commutators 

(1 + a) [d-A{x),i;+{y)] 
(1 + a) [a_A(x),V^Uy) 



o -I- 

—6(5 (x — y)'tp\_{x). 



(6.54) 

(6.55) 

(6.56) 
(6.57) 



These properties show that the A field can be taken for the specification of physical states via the 
condition 

{phys'\A{x)\phys) = (6.58) 



just like in the ET formalism [35|, |^6| . However in the perturbation calculations, the property 
of free propagation is far less important than the presence of q-number commutators which can be 
nontrivial obstacles for Wick's contractions. Therefore when one works with A, then one should take 
other fermion fields 



(6.59) 
(6.60) 



which already commute with A. However, for these new field operators the Hamiltonian density 
operator changes drastically (the field (j) decouples from the fermion currents) 

2 



ntZi = I (d.C,)'' - ad,Kd,4> + J faA + ^ * J+] - d,C,^ * J+ + ^S^IcT^ * ^lc, (6.61) 



d- 



d- 



^^^^ id 
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where now we have 



J+ = ^/2xh, (6.62) 
Clc = [-idia' + M(] + eaa')x , (6.63) 
eic = {id^x^a' + Mx^p)+ex^a'Ci. (6.64) 

One can verify that this new system describes the QED for the LC-gauge condition and the pertur- 
bative gauge field propagator would have the CPV-prescription for spurious poles. On the one hand 
this shows how different gauges can be linked together at the level of quantum field operators, while 
on the other hand this explains how easily the causal poles may be replaced by the CPV ones. 



6.2.1 Perturbation theory 



The perturbation theory is formulated in the interaction representation with the free Hamiltonian 
density 



Tin 



o^+^ = limWf+^ = ^X' + l{d,af-ad.4>dX+^d:;^H^ 



where 



and the interaction Hamiltonian density can be divided into two parts 

o/D+l _ -riD+l _ ijD+l _ njl I -1/2 



where 



Hi 



1 



int 



1 



1 



^2 " \id- - ed-(j) id- 
e 



+ -^tpla'iCi + dicP)— ^r- * Co + -^^0^7^ ^^ 

a/2 id- — eo-cp V2 lO- — eo- 



* (Q + ^^<P)^^J^ 



=2 

V2' 



1 



n 



int 



dia + 



a 



l + a 



-A 



id- — ed- 
a 



:a' * id + di<P)^l^^ 



l + a J d- 



A ) -1 * J+ + 1 



2 \l + aj \2d 



f 1 



(6.65) 

(6.66) 
(6.67) 

(6.68) 



(6.69) 
(6.70) 



We notice a close analogy between the present case and the model from Section |2[ therefore here 
we will point out only the differences between these two cases. The first one comes from the free 
equations of motion 



[{l + a)d+d--a^i_]\ 
[{l + a)d+d- - qAx]0 



(6.71) 
(6.72) 



where here we have the multipole non-covariant dynamical field (p which has no true Fourier repre- 
sentation. However, from the commutation relations ( |6.46 ), ( |6.47 ) and (6.48) one can derive the form 
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of chronological products 



(0|TQ(x)C,-(y)|0) 
(O|r0(x)A(y)|O) 
{0\T(^{x)(^{y)\0) 

Ti^+{x)ijl{y) 



SijDpix - y) , 
-Gipix-y) , 
a^Glp{x-y) , 
iV2A+dlAF{x-y,M^) 



(6.73) 
(6.74) 
(6.75) 
(6.76) 



where the covariant Dp{x) and the non-covariant G^p{x) massless Feynman propagator functions are 
defined in Appendix |B.1| . The combinations of independent modes which are coupled with fermion 
currents are 



A, 



1 

— * 

Gi + d, 



diG,. 



-A + ^ A I 



1 + a 



1 + a 



(6.77) 

(6.78) 
(6.79) 



with the chronological products given by 



{0\TAi{x)Aj{y) 
{0\TA_{x)Aj{y) 
(0|ri„(x)i_(y) 
{0\TA+{x)A^{y) 

{0\TA+ix)Ai{y) 
(0|ri+(x)A+(y) 



0) 
0) 
0) 
0) 

0) 
0) 



6ijDp{x -y)+ a'dfd]Glp{x - y) , 
a^dldyGlp{x-y) , 



a 



-^dldlGlp{x-y) 



(6.80) 
(6.81) 
(6.82) 

-I , -arv" ■ ^-ai-v- ^, ; (6.83) 

1 + a L J 
i(l - a)dldy{Gip * Dp){x -y)+ a^dldfGlp{x - y) , (6.84) 



a 



Gip{x-y) + a^A^Glpix-y) 



i2{l - a)dldl{Gip * Dp){x -y)+ a^dldlGlp{x - y) 



+ 



1 



(1 + a)2 dl 



* 5(x - y) , 



(6.85) 



where no explicit inverse Laplace operator appear; thus the limit d — > 2 can be taken already at the 
level of independent modes and the free propagators in 3+1 dimensions have the Fourier representation 



{{)\TA^{x)A,{y)\{)) =i 



(27r)4 



-ik-{x—y) 



1 



2k+k. 



— k± + ie 
c^ky k^ 



-9t^u + 



fc+(l-a)[<^fc^ + n;;^fc^ 
(1 + a)k+k- — ak\ + ie' 

1 1 



[{1 + a)k+k^ -akl + ief (1 + a)^^^"^ kl \ ' ^^'^'^^ 



where n^*-" = l,n^ = n^^ = 0. However, the Wick contractions of gauge fields and direct interac- 
tions of currents will give rise to the modification of perturbative propagators 



- y) = ^ (0 \TA,{x)A,{y)\ 0) + j^^^ * - v) 



(6.87) 



which no longer have any spurious pole with CPV prescription 



d^k 



-ik-x 



2k+k. - A;2 + 



-9iJ.u + 



k+{l - a)[n'^^k^ + n'^^ky 
(1 + a)k^k^ — ak'j_ + ie' 



oPkyk^ 



[(1 + a)k+k^ - ak^ + ie]"^ j ' 



(6.88) 
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The cancellation of non-covariant terms during tfie Wick contractions of fermion fields is identical with 
that in Section therefore, quoting those previous results, we say that the perturbative calculations 
based on the interaction Hamiltonians ( |6.69 ), ( |6.70| ) and canonical free propagators are equivalent to 
those with covariant interaction vertices and covariant perturbative propagators. In the perturbative 
gauge field propagator ( 6.88| ) with all causal poles we can take the appropriate limits to the Lorentz 
gauge (a 1) 



and to the LC-gauge (a 0)^ 



-ik-x 



1 



2k+k^ - kj +ie 



-9fiu + 



kykit 



[2k+k^ - akj_ + ie] 



(6.89) 



—ik-x 



2k+k- -ki +ie 



k+[nl;^k^ + n]fku\ 
k+k^ + ie' 



(6.90) 



When we compare ( |6.89| ) with the double limit (a — > and m? — > 0) of ( |2.96| ) we find that they 
are the same and also coincide with the ET propagator. The expression ( |6.90| ) represents the ML- 
prescription for the LC-gauge within the LF procedure. Contrary to other attempts |37], |38] which 
used two front surfaces for free gauge field sector, the present results are valid for the perturbative 
QED with fermions. 



7 Electrodynamics of charged scalar fields 

In previous sections we have discussed electrodynamics of charged fermions which describes the phe- 
nomena where electrons interact with photons. However it is also interesting to consider charged 
matter consisting of other fields. In the usual ET approach, scalar fields are treated as the simplest 
possible choice for matter fields. In the LF formalism, electrodynamics of scalar fields is far from 
being trivial. Below we will discuss only one choice of gauge fixing condition, the LF Weyl gauge 
in 1+1 dimensions, which has particularly strange properties. Electromagnetic currents built from 
scalar fields contain derivatives of matter fields j^"-'-"-^ = ecf)^ dfj. 0, which are dramatically different 
from the fermion currents. At LF this means that the complete dynamics of scalar fields directly 
manifests itself in coupling with the gauge fields. While such phenomenon has been recognized a long 
ago m, it is only here that a solution to this problem has been found. 



7.1 LF Weyl gauge in 1+1 dimensions 

When the LF Weyl gauge ^+ = is strongly imposed on gauge fields, the electrodynamics of scalar 
fields is described by the Lagrangian density 

^wlyr'"" = \ {d+A.f + {d. - ieAJ) 0ta+</. + 9+0t (a_ + i^AJ) <t> - m^cpU (7.1) 

and its equations of motion (Euler-Lagrange equations) are 

dlA^ = -ie(^(p'^d+(l}-d+(pU) , (7.2) 

2(9_ +ze^_)a+</> = -ie09+^_-mV, (7.3) 

2 (a_ - ieA_) d+^^ = iecjy^d+A^ - m'^cj)^ (7.4) 



^The ML-prescription is written here in the equivalent (in the sense of distributions) form 
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This means that aU fields are dynamical variables, but as usually, the Gauss law 

G = d-d+A^ - ie (cji'^d-ct) - </.5_</>"f) + 2eVV^- (7.5) 

is missing and has to be added as an extra postulate of the physical quantum theory. The canonical 
momenta conjugated to all dynamical fields 

= d+A. , (7.6) 
= (5_ -ie^_),/.t , (7.7) 
n^t = (a_ +ie^_)</) , (7.8) 

show the presence of the second class primary constraints (according to Dirac's nomenclature) which 
are characteristic for LF dynamics of relativistic particles. However, contrary to our previous treat- 
ment of similar constraints, due to the presence of other fields we cannot discard (7.7 and 7.8) as 



superficial constraints. This makes a sharp distinction between these expressions and canonical mo- 
menta for Ai in QED or even for free scalar fields. In our case, Dirac's method of quantization for 
constrained systems would give nonzero brackets (commutators) for scalar fields and 11" or between 
two momenta 11". All these brackets and also the one for scalar fields would depend functionally on 
A-. Thus everyone should understand those who decided to choose the LC-gauge A^ = where all 
above problems totally disappear. 

Here we follow our method of dealing with constrained system where we pay special attention to 
equations of motion. Eqs. (7^7.4) can be rewritten as the first order equations in 5+ 



d+ir = -ie[fd+(t)-d+(t)^(t)j , (7.9) 

2 {id- - eAJ) d+4> = e(t)ir - im^cl) , (7.10) 

2 {id_ + eA-) d+(j)'< = -e(l)^U- - im'^4>^ , (7.11) 

= H-. (7.12) 

Then using the Green function {id- + eA-)^^[x^ ,y^]'^^ for the covariant partial derivative 

{id- - eA.f {id. - eA_)-^[x- ,y-] = - {id. + eA.y {id- - eA.)-^[x- ,y-] =6{x~-y-) , (7.13) 

we can transform these equations into the form of Hamilton equations of motion 

d+cl){x) = ^Jdy-{id.-eA-)-^[x-,y-][ect)U--im^^){x+,y-), (7.14) 

d+(t>^{x) = ^ J dy' (^ecp^W +im^(t)^^ {x+,y-){id- - eA-)-'^[y-,x-] , (7.15) 

d+Il~{x) = ^(p^x) J dy' {id- — eA-)^^[x^ ,y^] (^—ie(l)Il~ — rn^cj)^ {x^ ,y~) 

+ |(/>(x) J dy^ (ie0l'n~ - m^^"^) (x+,y")(i5_ - eA-)-^[y- , x'] , (7.16) 

d+A-{x) = U-{x) , (7.17) 

or in the self-explanatory matrix notation (which will be used hereafter) 

d+(l) = ^{id- -eA-)-^ * (e(j)U- -im^(l?j , (7.18) 



26 



The perturbative definition of this Green function and the n otat ion for convolutions of integral operators, which is 



specially suitable for the further analysis, are given in Appendix 



B.2. 
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lie, 



^(t)\id- - eA-y^ * (-fe^n^ - + | (ie(l3'^ir - m^^) * {id- 



ir. 



eA^ 



(7.19) 

(7.20) 
(7.21) 



Now we encounter a real problem because, on the one hand, these equations evidently describe 
interacting fields, while on the other hand, the canonical Hamiltonian density has a free field form 



njQED 



U-d+A- + U^^d+cP'^ + U^d+<)) - Cweyi = :^{^'? + m 



(7.22) 



and does not depend on the coupling constant e. Even more, the generator of translations in the 
direction x~ is not a kinematical operator but depends on interactions 



dy- \ird-A- + (9_ - ieA-)(j)d-(y + (5_ - ieA-)(j)'^ d-(i) 



(7.23) 



These observations are apparently in a conflict with the original Dirac analysis of the front form of 
dynamics . The only way out of these problems is to allow the Dirac brackets to contain interactions 
and it really happens here because we have non-vanishing expressions 



{n {x+,x ),A_(x+,y )}^5 
{ n" (x+ , X" ) , (/)(x+ , y " ) } 

{n-(x+,x-),(/.t(x+,y-)}^^ 

1 0(x+ , ) , (/)"^ (x+ , ) } 
{n-(x+,x-),n-(x+,y-)}^^ 



- 5{x 

e , , 



y 



x^ ,x ){id^ 



eA. 



[X ,y J , 



(p^{x~^,x ){id- — eA-) '^[x ,y ], 
{id- - eA-y^[x' ,y'] , 



-i — (p\x~^,x ){id- — eA-) ^[x ,y ](l){x~^,y ) 



+ i^(f)''{x^,y ){id- — eA- 



[y ,x 



(X ,x 



(7.24) 
(7.25) 

(7.26) 
(7.27) 

(7.28) 



Now it is not difficult to show that these brackets lead to the correct equations of motion for all field 
variables ( [7.16 -7.15) and also give the expected translations in the direction x~ 



{n-(x),p_(x+)}^^ = d-u-{x) 

{A-{x),P-{x+)}^^ = d-A-{x) 

{0(x),P_(x+)}^^ = 5_0(x), 

{0t(x),P_(x+)}^^ = 9_0t(x). 



(7.29) 
(7.30) 
(7.31) 

(7.32) 



If one would like to take these brackets as a basis for respective quantum commutators, while defining 
the canonical quantum theory, then one would end up with a hopeless problem of the perturbative 
calculations. When the interaction is not located in the Hamiltonian but in commutators, then the 
definition of the interaction representation is not unique, if at all possible. 

7.2 Dressed scalar fields 

In our previous analysis we have found that LF dynamics allows for using various forms of fields which 
differ in both the evolution equations and the commutator relations (e.g. the Lagrange multiplier 
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fields in Section |6|). Thus we expect that there are also various forms of scalar fields, one of them 
may have free commutators. Usually at LF, the dressed scalar fields have the field-dependent phase 
factor like 

^= exp|-i e(a_)"^ *yl_|(/) , (7.33) 
= ^texp {i eA- * {d-Y^^ . (7.34) 

However this really means a gauge transformation to the LC-gauge and evidently is not a true solution 
for the LF Weyl gauge problem. Therefore we take another possibility and define a new scalar field (/? 

= >V_i[a]*(^ , (7.35) 

where the integral operator W_i[a] is defined in Appendix |B.2| . For simplicity, the other scalar field 
(j)^ is not changed]^ Now one can easily check that for the pair of scalar fields ip and their Dirac 
bracket has already a free form 

{v?(x+,x-),</<t(x+,y-)}^^ = -'-{id_)-\x- -y-). (7.36) 
Encouraged by this result we find the scalar field contribution to the Lagrangian 

j^scalar ^ q^^] ^ ^ ) + (9_ _ ^e^_ ) O+cf) - cf) = 

= d-cj)^ * dj^Lp + d^cj)^ * d^Lp — m^cj)^ * W-i[a] * ip 

+ie(j)^ * 5+A_W_i[a] * p (7.37) 

and then calculate the canonical momenta conjugated to scalar fields 

= , (7.38) 

n^t = d.p, (7.39) 

which are already free. This means that we have again trivial primary constraints for these momenta 
which can be discarded in the Hamiltonian approach. Also the third momentum (for gauge field) 
has been changed 

n = d+A^ - d-A+ + ie(l)^W-i[a] * p , (7.40) 

but there are no constraints connected with this mode, so this change introduces no problems into 
further analysis. Now the complete Hamiltonian for the LF Weyl gauge is 

H?an" = ^ J dx- (U-ie^^W^i[a]*py +m^^UW-i[a]*p (7.41) 

and evidently it contains interactions, while the translation generator 

P_ = iy"(ix~H*9_^„ +2a_0U5_(/9 (7.42) 

is kinematical. Finally, we give nonzero commutators at LF 

2dl\(j)\x+,x-),p{x+,y-)] = i6{x- -y-) , (7.43) 
[U{x+,x-),A_{x+,y-)] = -i6{x--y-), (7.44) 

which also have canonical free forms. Therefore we see that the redefinition ( [7.35| ) of the scalar field 
leads to the canonical description of scalar QED at the LF. 

^'^Also another choices of dressed scalar fields are possible, specially the symmetrical choice for both fields. These 
issues will be discussed elsewhere. 
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7.3 Perturbation theory 



Now we can define the interaction representation with the free evolution of field operators which is 
given by the free Hamiltonian Hq 



Ha 



dx n + m 



(7.45) 



The evolution of states will be given by the interaction Hamiltonian Hj 



+ m^(/)J * (W_i[a] - 1) * 



(7.46) 



At low orders (in e) of the LF perturbation, different contributions generated by the above Hamil- 
tonian formally sum up to the covariant results, which can be obtained within the ET approach]^ 
However, if we want to prove the equivalence of both perturbations for all orders then we need to 
follow the methods used earlier in Section |2[ In order to separate the effects of scalar and gauge field 
contractions, we divide the interaction Hamiltonian into two parts 



Hi = Hj+Hj, 

H^i = \j dx- (fe4w_i[a] * = j dx^ J^. 



The functional form of the Wick theorem 



Texp —i J aJn = : exp —i J Jn + • ~ ^xp — - J : exp —i J aJjj 



(7.47) 
(7.48) 

(7.49) 



(7.50) 



indicates that the Hamiltonian Hj can be equivalently substituted by the linear term / crJn, where 
the new field a has non-zero Wick's contractions 



(0 \T+a{x)a{y)\ 0) = i6'^{x - y). 
Then from the canonical propagators 



(o|r+n(x)n(y)|o) = 

(0|T+^_(a;)n(y)|0) = 
we find the Wick contractions for the modified field H" 

;o|r+n-(a;)n-(y)|o; 
:o T+A^{x)ir{y) o; 



E\{x - y) , 

= n- -Fcr 
i5^{x - y) , 
E^x-y). 



(7.51) 



(7.52) 
(7.53) 



(7.54) 
(7.55) 



Thus effectively we can take the equivalent interaction Hamiltonian which is bilinear in scalar fields 



H'^f^ = (/)q * |— ieH -I- m^aj * yV_i[a] * = cj)]^* Hcov * VV_i[a] * (f. 



Having the Wick contraction for the pair of scalar fields 



T 



(p^ {x)(p{y) = Af{x - y,m'^) 



(7.56) 



(7.57) 



^*The ET interaction Hamiltonian also contains one noncovariant term which cancels the noncovariant part of the 
second time derivative of gauge field propagator [p9[. 
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we can write the functional form of the Wick contractions for scalar fields 
r+ exp -i * Hcov * W_i [a] * c^j = 

= exp ~* I ("^^ + ^ * * -^coD * W-i [a] * + AiT * ^ 
= : exp —i l^(f>^ * Hcov * W'-i[a] * (l + iAp * Hcov * W'_i[a]^ * (^j 

X exp-Trln(l + iAi;'*#co,;* • (7.58) 

In Hcov there is the integral operator m?{id^)^^ which has a convolution either with the field (j^ 

4 * rr?{id-)-^ = -2id+(pl (7.59) 

or with the propagator function Ap 

Af* m^{id-)-^ = -2id+AF d+ +i{id-)-^. (7.60) 
Now it is quite easy to check the following factorization: 

l-iAp* Hcov * W-i[a] = (1 - iAp * Hcov) * W_i[a] , (7.61) 

where 

Hcov = -ieU - 2ie d+ A_ , (7.62) 
and next we can write the expression for the effective scalar field contractions 

exp —i ^(j)^ * Hcov * VV_i[a] * (/?| = : exp —i ^4>^ * Hcov * (1 + Ai? * Hcov*y^ * ^jj x 

X exp -Tr In {I + iAp * Hcov) x exp -Tr In W_i [a] . (7.63) 

In this way, we have found the effective Feynman rules for scalar and gauge fields in the usual form 
of perturbative propagators (|7.57|) and ( ^.511) , respectively. The vertices are given by 

Hcov = -ied+A_ - 2ie d+ A„ (7.64) 

so at any vertex we have the factor (2p+ — fc^), where A;+ is the momentum of gauge field line and 
P-\- is the momentum of scalar field line at its <f)^ end.0 It is very interesting that here for the scalar 
field QED all noncovariant contributions boil down to the closed loop factor 

exp — Tr In W-i [a] ~ exp Tr In [id- — eAJ\ (7.65) 

and according to the same arguments that we have used for the fermion field case, we can omit these 
contributions completely. 



^^Due to the momentum conservation at each vertex, this is equivalent to the symmetrical rules where momenta of 
both scalar lines are taken, for example see |45[. 
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Part IV 

Finite volume QED 

8 LF Weyl gauge 
8.1 The DLCQ method 

Notation for models in a finite volume of LF are generally given in [40|. Below, for convenience, we 
present some of the basic steps of the DLCQ method underling these points which are different from 



We choose the space of LF as a " hypertorus" : —L < < L i —L± < x± < L± and impose periodic 
boundary conditions for boson fields. For fermions we choose antiperiodic boundary conditions and 
this choice means that, contrary to bosons, fermions will solely have modes with all nonzero compo- 
nents of momentum. 

Zero modes of boson fields can be discussed using the classification introduced in |41]. We denote the 



full gauge field as V^(x) and this can be expanded in Fourier modes. We distinguish, respectively, the 
simple zero mode 

(8.1) 



-L 2L 



and the normal mode gauge field: 



The latter degrees of freedom are known to represent the usual propagating photons in the light-cone 
representation, and as such we reserve the symbol to denote them. From the simple zero mode 
one can build the totally space-independent global zero modes 

f^^ dx-£x^ 



L 



where in future we shall write d'^x for dx~d'^x± and suppress the limits of integration. Evidently, the 
are 0-1-1 dimensional fields, namely quantum mechanical variables - thus the notation q. Finally, 
the simple and global zero modes can be used to build modes with no x~-dependence but no constant 
part in x±, i.e. the proper zero modes: 

r+L fj^rj.- 

af,{x±)=l ^ -^V^{x- ,x±) - q^ . (8.4) 



-L 



When the decomposition of modes is complete we shall refer to one of the normal yl^, proper zero 
mode a^ and global zero mode sectors. 

For all the above sectors one needs the corresponding delta functions. We adopt the notation that 
the periodic three-dimensional delta function is represented as 6^^^ {x — y), which includes the zero 
modes. For the antiperiodic delta function, a subscript 'a' is appended: 6a (x — y). The explicit 
difference between these two objects can be easily seen by expanding in discrete Fourier modes. Next 
we distinguish the delta functions appropriate for each mode sector for periodic functions. Thus in 
the normal mode sector we must subtract the x~ independent part of 5^'^^ , and so define 

5(3) (f _ y) ^ 5(3) _ y.) _ (^^ _ . (8.5) 
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In the proper zero mode sector we must subtract the overall two-dimensional volume factor in defining 
the relevant delta distribution: 



x± - y± 



1 



(8.6) 



When reaching a canonical formulation, the QED Lagrangian, expressed in terms of the complete 
fields V^, if: and ij^^ takes the standard form 



L 



(fx 



^7) 



Once the boson field is decomposed into its different sectors, V^(x) = ^/^(x) + a^{x±_) + (7^ , the 
Lagrangian Eq. (|8.7D breaks into three parts 



nm ~ ^pzm ~ ^gzm i 



where 











Lpzm — 


J (fx 


-\{^^lau - 






Lgzm — 


{8LLl 




- eq^ / d^xi^-f^'ip . 





(8.9) 
(8.10) 
(8.11) 



These formulas in turn can be simplified by decomposing the total electromagnetic current, = 
—eip'j'^tl^, into its normal mode and proper and global zero mode parts: J''^ = + J^zm + 
Fermion current, while being bilinear in fermion fields, satisfies periodic boundary conditions and its 
global zero modes are given by 



Q' - 

The proper zero modes are 



3^ 



8L L\ 



^'±{x)i\)±{x) 



dh 



8L Li 



V'+(x)a*V'-(^) + V'l(^)a*V'+(^) 



iV2 



dx 



'^pzm (-^-L ) 



L 2L 
^ dx 



■vl(x)v±(f)-g 



± 



2L 



(x)aV- (x) + V"! (x)q V+ (x) - Q 



t (;^\^i. 



Finally we can give the normal modes for fermion currents as follows: 
Jnmix) = -e\/2V'i(f)^±(f) - J^2^(x_l) - , 



ijl{x)a'tlj-{x) + ij'L{x)a'ij+{x) - J' {x±)-Q 



(8.12) 
(8.13) 

(8.14) 
(8.15) 



(8.16) 
(8.17) 



This notation is diflterent from m» in order to avoid collision with external currents j'^ and fermion currents J'^ 
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8.2 Canonical formalism for LF Weyl gauge 

The Weyl gauge condition can be imposed strongly, namely = a+ = g+ = at the classical level. 
As usual, this means that Gauss' law (given explicitly later) appears as a constraint to be imposed on 
the physical states. Our procedure for carrying the canonical procedure is as follows: we will analyse 
different subsystems, where only one field sector is treated in terms of independent degrees of freedom 
while the remaining fields are regarded as non-dynamical external fields and/or currents. Then we 
will exchange non-dynamical modes for effective interactions of dynamical ones and will give a simpler 
(though nonlocal) Lagrangian where only dynamical fields are present. This will result in a sequence 
of equivalent effective Lagrangians which contain fewer modes but have the same Euler-Lagrange 
equations as those which are generated by the primary Lagrangian, provided the constraint equations 
are implemented for non-dynamical fields. This procedure is based on the observation that different 
Lagrangians can lead to the same system of Euler-Lagrange equations though they may have very 
different constraint structure. One can feel free to choose the most suitable one for carrying out the 
canonical quantization procedure. 

8.2.1 Proper zero mode sector 

The Lagrangian Eq.( 8.1C| ) can be written explicitly in LF coordinates 



1 1 



(8.18) 



where the Weyl gauge condition a+ = has been explicitly imposed. Because we are interested 
here in boson fields we will treat fermion currents as arbitrary external currents, however without 
introducing any distinction in notation we hope that this will cause no misunderstandings. Thus the 
Lagrangian ( ^.18 ) leads to the classical equations of motion 



d\a- = did+ai + Jp^^ , (8.19) 

-did+a^ = A±ai - didtak + Jp^rn > (8-20) 

which impose the parameterization of the field Oj 

^ * (jLm + i^ivr ) . (8.21) 



A, y P'"" ' 2L 
This new field tt has the Dirac bracket 

{7r(xx),a_(yx)}DB = -Sl?Hx± - y±) , (8.22) 

which is the only nonzero bracket in this sector. The canonical Hamiltonian contains the effective 
nonlocal terms 

2L fir 1 ■ f 

^wTyl = '^^^^■^pzm^*Jpzm- J C^^^^^TT— * J^^„ - 2L J (fx^a^J'^^. (8.23) 

One can check that the effective equations of motion which follow from the above Hamiltonian and 
bracket agree with the Euler-Lagrange Eqs. (|8l9| , |]2|). Therefore our Hamiltonian and brackets 



describe the same classical system as the primary Lagrangian Eq.( 8.1q ) and one can give an equivalent 
Lagrangian density 

^eff = ^^^+^- ~ ^d"" = ^^^+"- + 2L^A7 * ^^'^P^™ + O'-'^pzm - 2'^pzm-^ * Jpzm > (8-24) 

which directly leads to the correct bracket and equations of motion. 
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8.2.2 Normal mode sector 

In the second step, we analyse the sector of normal modes of gauge field potentials treating normal 
modes of electromagnetic currents as arbitrary external sources. From the Lagrangian ( |8.9D we 
take these terms which contain 

£^rey/ = 9+ A, {d^A, - diA^) + i {d+A^f - ^ {diA, - d^A^f + A.J- + A,j' , (8.25) 
and find the Euler-Lagrange equations of motion 

d+{d+A.-diA,) = J-^, (8.26) 

{2d+d--A^)Ai = d^{d+A^-djAj) + r^:^. (8.27) 

We see that these equations have the same structure as the respective equations in Section ^ in the limit 
a — > 0, therefore we can adopt our previous results here. Thus we take a new field 11 = — diAi 

and get the diagonal structure of Dirac brackets 

{A.{x),U{y)}DB = S^nHx-y), (8.28) 
2d^.{A{x),A,{y)}^^ = -6,j6l^){x-y) , (8.29) 

while all other brackets vanish. The canonical Hamiltonian 

n^eyi = + mA, + \ {d^A^f - A^j;;^ - yl, j;^ (8.30) 

generates equations of motion which are equivalent to the previous ones and can be used for defining 
an effective Lagrangian 



£^,yi = d+A;d^A, + Iid+A.-n^,yi 



d+A,d-Ai - ^ {diA.f - -{lif + n {d+A^ - d^Ai) + a.j-^ + Aai^. (8.31) 



Having analysed the canonical structure of the gauge field sector, one can substitute the Lagrangian 
Eq.( ^.ldD by Eq.( |8.24| ) and the boson part of Eq.( |8.9D by Eq.( ^.31 ) and instead of the total Lagrangian 
Eq.( p.8p , one can work with the effective Lagrangian 

C^iyi = d+Ad.A-^{diA,f -^{U)' + U{d+A.-d,Ai) + ^7rd+a- 

- MV'+"^7V- - MV'-"^7V+ - eV2iljJ^.V. (8.32) 

- eV{ (V'+'^aV- + '0-"^aV+) - \jpzm * • 

In this expression, 

= Ai-^d^^*7r + qi, (8.33) 



namely it is the original Vi but with its proper zero mode Oj expressed as in Eq.(8.21) and the new 
"current-current" interaction subtracted]^ The decomposition of V- remains unchanged. 

^'^This term has been explicitly reintroduced as the last term in Eq.(B.32). 
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8.2.3 Fermion sector 

In the next step, we take the fermion part of the effective Lagrangian Eq.( 8.32] ) 



'-Weyl 



(8.34) 



The non-dynamical modes are the fermion field ip- and the global zero mode Qi. These are determined 
by the following differential equation and global integral condition: 



= Q' 



8LL^ 



The first equation leads to 
1 1 



qi)a 



id- — eV- 



where 



* vr (x) 



ij+{x) 



(8.35) 
(8.36) 

(8.37) 
(8.38) 



Note that these are not yet the solutions for the dependent fermion field ip- because these fields 
appear also on the right-hand side in the zero mode currents Jpzm- However, one can introduce them 
into the definition of J^^„ + Q' (see Eqs.( |8T^j8l5D ) 



2L 



2L 



* Jpzm ) {X±) - Qk 



where 



r(x^) 



1 

V2 
1 



dx ilj\_{x)a 



1 



-|- I dx ^\x)a^ 



id- — eV- 
1 

id- — eVL 



* -0+ I (x) , 



M'^ix±) = 5'^V2 J dx^il^lix) (^ .^ * (x) = 6'^M'^ix^) . 



Now, from the constraint = one gets the differential equation 



A, 



^M\x,) 



which has a formal solution 



ix±) 



e 

'2L 



(8.39) 



^.40) 
^.41) 



.42) 



^.43) 
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in terms of the functional G(±)[x±,y±; Ai'^] introduced in Eq.(B.40). Finally, we may express the 
non-dynamical fermion field as 



^2 \id--eV. 
whereby we obtain 



* ^) + ^ {G(±)[M'] * r^) (xx) [ ^Q_^^y_ * (^) , (8.44) 



1 



id^ — eV- 



(x) + cfx^Vix^) (g^^)[M^] * r) {x^) (8.45) 



as the Dirac Hamiltonian for unconstrained fields. Just as before we can give the effective classical 
Lagrangian for the fermions 



id- — eV- 



^.46) 



8.3 Quantum theory 



Having eliminated all non-dynamical fields we may now proceed by substituting the fermion part of 
Eg. ( ^.32 ) by that given in Eg. ( 8. 46 ). We thereby obtain the total effective Lagrangian 



o2 



1 



1 



2L ' 2 

2l 



^.47) 



Jd- — eV- 

as the starting point for the canonical guantisation procedure. It generates the Euler-Lagrange 
eguations of motion which agree with the dynamical eguations of the primary Lagrangian Eg. (8. 7) 
with all non-dynamical eguations (formally) implemented. The canonical guantisation is simple here 
and one gets the egual-x^" guantum commutation relations 



[Uix),A. 



-i6l^\x-y) , 



{iji{x),^+{y)} = i= A+<5i3)(x-y) , 



[■^ix±),a-{y±)] 
2dl[A.,ix),A,{y)] 



V2 

X. . 



p yX± - y±) , 
i6ijd';^\x - y) , 



(8.48) 
(8.49) 

(8.50) 

(8.51) 
(8.52) 



and the guantum Hamiltonian, like Ceff, which comes from H^"'", ni"" and Hj^ 

1 



fer 



(fx 



^{U{x)f + mAi{x) + ^ {diA,{x)f 



+ 



-\2 



-ip-) 



+ 



1 



V2 



d^xi\x) 



1 



id- - eV- 



* e (x) + 



1 

22I 



IQLL\ 

dVr^(^±)(a(±)[A^Vr) (xx), 



^.53) 



tacitly defining the ordering. One can show that, due to the effective eguations of motion, the Gauss 
law operator 



G{x) = 9_H(x) + 2d-diAi{x) - A_l^-(x) - A_La-(x_L) - e\/2V^Y(f)V'+(x) 



.54) 
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is ^'''-independent. It leads to a classical first class constraint, G ~ 0; namely it must annihilate 
physical states in the quantum theory. Furthermore, it is intimately connected to the residual gauge 
symmetry 



A_{xf 

-h 



nhA^{x)n{ = Ai{x) - dih{x) , 
nhA_{x)nl = Ai{x) - d^h{x) 

Q.hii{x)n\ = n(x) + A_L/i(f) , 



Q.ha^{xi_)Q\ = a-{x±) , 

0.hTT{x±)^{ = ■k{x±) + A±_ J ^dy~ h{y',x±] 

fihP'^i =p~ , 



Jeh{x) 



where 



(8.55) 
(8.56) 
(8.57) 
(8.58) 
(8.59) 

(8.60) 
(8.61) 

(8.62) 



8.3.1 Translation Generators 



With the effective Lagrangian density Eg. ( ^.47 ) one can calculate the canonical momentum-energy 
tensor 



rpUfl 



+ 



S {d.Aij 

1 6Leff 



'6 {d^i^+) 2L5 (d^TT) 



2L5{dya_) 8LLi5{d^q. 
Then from the generators of translations = J d^x T^^{x), the spatial translations are 

= j -d-AkdiAk - lidiA. - ^^/2^/;t9^V'+ - ^^d, 
d^x \d-Akd^Ak + lid-A. + iV2^l)\d-i}^ 



^.63) 

^.64) 
^.65) 



Now from the (anti) commutation relations Eqs.( ^.48 -g.50) one can recover the correct Heisenberg 
relations for all dynamical quantum fields ipj = (^fc, H, vr, a_) 



^.66) 
^.67) 



and this confirms the translation invariance of QED in the Weyl gauge. We note that the generators 
and are not invariant under the residual gauge transformation with gauge function h{x): 



Pi 



pi 



P^ + J d^x G{x)d-h{x) 
d^x G{x)dih{x) , 



P 



(8.68) 
(8.69) 



and this is connected with the lack of gauge invariance of the canonical energy- momentum tensor. 
We return to this below. 
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8.3.2 Implementation of Gauss' Law 



As we have mentioned earlier, the physical quantum gauge system has to satisfy Gauss' law |42]. 
However this cannot be implemented strongly as the condition for quantum field operators which 
would be incompatible with the commutation relations but rather as the condition for physical states 

G(x)lphys) = 0. (8.70) 



In 1 40 1 the method of quantum mechanical gauge fixing |^3[ has been used and details of this formalism 
are given therein. Here we will give only some most important results. Two gauge transformations 
are defined with the help of the Gauss law operator 

Ui[^] = exp(^-i J d^xg{x)^[x;A_]^ , (8.71) 

= ex.p (^i J cPx±p2{x±)r][x±;iT]^ , (8.72) 



where 



^[x;A.] = (^J-*A_)(x), (8.73) 
g{x) = G{x) - a_n(x) = 2d-diAi - A^A^ - A±a- - e\/2'0+V'+ , (8.74) 

v[x±;tt] = (2^^) ' (8-75) 

= dx-iii^+ . (8.76) 

They allow us to express the Gauss law condition ( |8.7[1| ) in two pairs 

niphys') = 0, (8.77) 

nT 

— <7(f)|phys') = 0, (8.78) 
where 

Iphys') = C/i|phys), (8.79) 

and 

a_(xx)|phys") = 0, (8.80) 

giphys") = J (fxp2{x^)\phys") = , (8.81) 

where 

Iphys") = C/2|phys'). (8.82) 

Here we have the neutrality condition ( |8.81| ) as the only one constraint which survives from the infinite 

number of primary constraints which have no physical meaning. 

Also the physical Hamiltonian can be defined as it acts on physical states 

U2UiHUM\phys") = HY,:%hys") (8.83) 
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and due to ( ^.77 and |8.8[1| ) the field operators 11 and a_ can be omitted as cyclic variables 



H 



Weyl 
fin 



I ( djAj - V2e^ 



d. 



id- — e 



*x \ (x) 



+ 



2 2L 



.84) 



The operators T" and M" follow directly from the respective operators T and M in which ^ is replaced 
by 

X = [m7o - ia'di + ea'Ai]^+ (8.85) 

and V- by 

^ '3.86) 



9- - -r— * P2 



where 



^.87) 



This redefinition of q'_ can be viewed as an infinite renormalization of the field, which does not change 
the commutators. 

Similarly for the generator of other translations we have 



UiUiP'ulul 



- J d^x (d.Ajd^Aj + iV2i/jldii^+) = P},^ , 
U2UiP+ulul = J d^x {d^Ajd.Aj + iV2'il;ld-i;+) = p}i^ , 



i.89) 



which means their invariance in the subspace of physical states. 

Also we stress that in the physical operators, the dependence on the proper zero mode a_ has totally 
disappeared which contradicts the naive interpretation of this mode as a gauge-invariant physical field. 



9 Covariant gauge 

In the finite volume LF the Lorentz covariant gauge condition d^V^^ = would be equivalent to three 
independent conditions for each sector 

r d+A-+d-A+-d^Ai = {) 
d^V = ^ <^ d+a- - ditti = (9.1) 
[ d+q- = . 

However the last condition for global zero modes is not attainable by means of any gauge transfor- 
mation because mode q- is gauge-invariant]^ Therefore here as the Lorentz covariant gauge we take 
the following conditions: 

d+A-+d-A+-diAi = 0, (9.2) 

9+a_ — diQi = , (9.3) 

g+ = 0, (9.4) 

■^^This resembles the situation in the LC-gauge in the sector of zero modes, where conditions a_ = g_ = are not 
attainable as gauge conditions mm. 
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with the LF Weyl gauge for global zero modes. In Lagrangians for different gauge field modes the first 
two conditions will be implemented by means of the Lagrange multipliers and Ap^m, respectively, 
while the third one will be imposed explicitly. 



9.1 Normal mode sector 

In the sector of normal modes we start with the Lagrangian density 

= {d+Ai-diA+){d-A-diA^) + ^{d+A^-d^A+f-^{diAj-djAf 

+ A-J-^ + A+J+^ + AiJi^ + Anm{d+A- + d-A+-d±A±). (9.5) 

Next we introduce the canonical momentum conjugated with A_ 

n- = d+A_ - d-A+ + Anm , (9.6) 

and then modify the Lagrange multiplier field 

^nm = Anm — -^g- * Jnm ' (9-'^) 

in order to separate dynamical equations 

(2a+5_-Ajn- = 25_J-^ + 5,4„, (9.9) 
(2a+5_-AJA = 5,A„^ + j;„ + ^* , (9.10) 



from the constraints 



A- = -^*d-[ll- +diAi-2\nm] , (9.11) 



diAi + Xnm - n + * J+„ 



(9.12) 



Prom the canonical Hamiltonian density 

= {d+A_)U- + id+Ai)W-£ = ^(n--Xnm-^*JL)' + \{diAj-djAif 

+ [Km + ^ * JL) diAi - d. [n- + diA - 2Xnm] ^ * Jnm ' AiJ^ (9-13) 

and the dynamical equations of motion we find the Dirac brackets 

2d-_{Aix),U-iy)}^^ = df5\x-y), (9.14) 

2a^{n-(x),n-(y)}^5 = A^5\x-y), (9.15) 

2dl{Ai{x),A,{y)}j,j, = -Sij6^{x-y), (9.16) 

2dl{Xnm{x),A{y)}^s = -df5\x-y), (9.17) 

while all other brackets vanish. However, another choice of field variables 

(j) = -^*{U- -2Xnm + diAi) , (9.18) 

d = Ai-di^*{U + djAj) = Ai-di(l), (9.19) 
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leads to a simpler structure of Dirac brackets 



2dl{C,{x),C,{y)}j,j, = -5,,5'{x-y), (9.20) 

2d'^{cl){x),XnmmDB = ^\x-y), (9.21) 
and dynamical equations are separated 

[29+9- -AJ,/. = A„„-^*J+^, (9.22) 

[29+9_ - A J A„„, = A+^* J+„ + a_J-^ + a,j;„ , (9.23) 

(29+9_-A^)a = j;^ + 9,^*J+^. (9.24) 

Next from the Hamiltonian density 

1 

2 yniiL -^L-^L ' "nm J ' ^ 

+ (A+(/> + did) (^Xnm + ^ * JL^ + (5_ J-^ + di4jj - r^^Ci (9.25) 

and the diagonal Dirac brackets (|9.20| , |9.2l| ) we can construct the effective Lagrangian density for the 
sector of normal modes 

J^Tv'"" = d.ad+a-UdiCj f -2d+Xnmd-(l) + diXnM-UXnm-:^*J^ 



1 

2d- 



(29,Q + A+0) * - (9_ J-^ + diJl\ + j;^C,. (9.26) 



9.2 Proper zero mode sector 

This sector is specially interesting because both its Lagrangian density 

= 5ia_ {dia+ - d+ai) + ^ {d+a^f - ^ {Oia, - d^aif + a^J^,^ + Ap,™(9+a_ - S^a,) (9.27) 
and the Euler-Lagrange equations 

d\a^ = di{d+ai - dia+) + J~ra - d+Apzrn , (9.28) 

-d+dia_ = A±ai - didja^ + Jp^^ + diAp^^rn , (9.29) 

= -A+a_ + J+^, (9.30) 

(?+a_ = ditti , (9.31) 

apparently indicate dynamical modes. However, a closer inspection shows that all fields are non- 
dynamical and explicitly depend on arbitrary external currents Jp^j^ 

^*(jP- + a^Ap,J , (9.32) 



"'^ - Aj. V'^A^ 

Apzm = * ip+'^pzm + '^iJpzm) ■ (9.33) 

Therefore there is no canonical structure here and this resembles the massive QED in Section 
All we can do now is to find the effective Lagrangian which properly couples the dependent gauge 
potentials ( |9.32| ) with external currents. The best choice is 

fpzm,eff ^ 7" T+ 4- - -8 A ^— ^— T f9 34^1 

'~'cov "^pzm ^^"pzm ' 2 P-^™ ^fJ.-^'-pzrn pzm ^^'^pzm ^^'^pzmi y^.Kj^j 

where the non-dynamical field Apzm is left in order to avoid the term 9+ J^^. 
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9.3 Fermion field sector 



Gathering our effective Lagrangians (|936|) and dpI) we liave tlie total effective Lagrangian density 



for QED in tlie Lorentz covariant gauge 

cm = d-ad+C,-\{d^C,f-2d+\nrnd-^ + d.^nn.d^4>-\\lrn-\{l^*J^ 

- eVi {ll)+^a'i)^ + V'_1'qV+) - ^J'pzm^ * Jpzm + Jpzm^ * Jpzm , (9-35) 

where the gauge potential is given by 

Vi = Ci + di(f> + * Apzrn + Qi , (9.36) 

V_ = + , (9.37) 

* {IdiQ + A^(t> - Xnm) +d+^ 



V+ = — * i2d^a + A^^ - Xnm) + d+— * Ap,m . (9.38) 



Though the equations of motion for fermion fields lead effectively to the condition Apzm = we would 
like to remove this field before analysing fermion sector. ^'^ This can be achieved by the redefinition 
of fermion fields 

V'i = exp (^-^e^ * Apzm^ V'± , (9.39) 



± = V'i exp ( ie— * Ap,m ] , (9.40) 



which in (9.35) changes gauge potentials into V^°^P^ 



vr = a+ 34 + Qi , (9.41) 

V^"" = d-(t) + q- , (9.42) 
vr = ^ * ma + A^4> - Xnm) . (9.43) 

Thus the non-dynamical field Apzm appears only in the term ^Ap^^i its equation of motion is 
Apzm = 0. Therefore we can completely disregard this field in the further analysis . 
The fermion Lagrangian density 



2d. 



Its presence would generate the extra primary condition which contains fermion fields 

Sdl 



34 



Strictly speaking these gauge potentials do not satisfy the Lorentz covariant gauge condition but rather 

dn, X rcov f T 



However the direct interactions of currents compensates the above non-vanishing current term. 
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- eV'"\ (^;^aV- + C^aV;) - ^4-^ * J'l.m + ^'p-^ * J'tm , (9-44) 

where all currents are built from the new primed fermion fields, generates constraint equations for 
dependent fermion fields V'^ and global zero modes qi 

{^V2d. - eV2Vl°- - * J'+ „^ = -ia'd^i,'^ + m7V^ 

+ ea^ (vr - * J'U) , (9.45) 

d^x (V'+aV- +'0'-aV+)(a;) , (9.46) 

with a similar structure as their counterparts for the LF Weyl gauge. Therefore we can incorporate 
those respective results from Section |8| and write 

(9.47) 

^ * 4- - ^« = * (9.48) 

Then we take the expression for currents 

J'Um = -^r- - j^M'\^)[M''] * r'^ , (9.49) 

+ j dx~i'\x) * .^^ ]_ aV;(x) , (9.50) 
7W''(xx) = V2 j dx~^'\{x)-^-^-^*^'^{x) , (9.51) 



where 



and 



= a_,/. + (7--eV2^* J'+^ , (9.52) 
^' = [m7° - ia^9i + ea' {Q + 9^0)1 V'' . (9.53) 



Next we find the effective Hamiltonian for fermion sector 



2 

nm 



■ J d^x^T'\x^)g^^) [M'^] * T'\x^) + e^/2V^; ' (9.54) 



1 e2 



2 2L 
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and finally we write the expression for effective total Lagrangian density 

- \^T'\x^)G^^)[M"]*T'\x^) . (9.55) 
9.4 Quantum theory and physical states 

When all constraints are solved at the classical level, then the canonical quantization is straightforward 
- one writes the (anti)commutation relations 

2dl[C,{x),C,{y)] = -i6,,5^^\x - y) , (9.56) 

2dl[<P{x),\nm{y)] = i5^^\x-y), (9.57) 

[i;\{x),^\{y)] = l=K+5f\x-y), (9.58) 

= -i, (9.59) 
while other relations vanish, and takes the quantum Hamiltonian in its canonical form 



+ 1^1 d'x^r'\x^) (a(x)[-M'Vr") (9.60) 

The presence of extra nonphysical modes (f) and A„m indicates that the Hilbert space contains also 
non-physical states. In order to select the physical states we will follow the analysis presented in 
Section |^ and introduce first the quantum field Anm 



Anm = ^nm + ^* J'nm, (9-61) 

which is a free field even for the fully interacting system 

(2a+a_ - Ax)A„^ = (9.62) 

and has nonzero commutators 

2d^_[m,Km{y)] = i5'i\x-y), (9.63) 

2d'L[Km{x),^\{y)] = e^'+{x)5^^\x - y) , (9.64) 



29^ 



Km{x),^\{y)\ = -e^%{x)5^^\x - y) . (9.65) 



The q- number commutators ( 9.64 ) and ( 9.65| ) would be very cumbersome in further analysis, therefore 
we introduce the dressed physical fermion fields 

ipphysix) = expie(p{x)'ip+{x) , (9.66) 
'^Ihysi^) = V'+(^)exp-ie(/>(x) , (9.67) 
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which aheady commute with A„m at LF. Next the quantum Hamiltonian can be expressed in terms 
of these new fields 



H, 



Lor 
quan 



+ 



1 e2 



-{P~ 



+ 



V2 



1 



phys * ^phys I 



+ 2 2L / ^^^Physi^^) [Q{^)[-Mphys]*^phys) , 



(9.68) 



where 



phys [-^-L J 



1 



id- — eV: 
1 



phys * ^Phys{x) 



^phys ^ 
( 

phys 



id- — eV_ 
1 



^a'i^physix) 



'^h^ * "^physix) 



and 



"^phys — ^'^phys'^phys 



V: 



phys 



phys 



^phys 



TT- * did , 







mj — la Oi 



a + ea'd 



'4'phys ■ 



Now we are in the position to give the condition for physical states 

{phys' \ Anm{x)\phys) = 0. 



(9.69) 
(9.70) 

(9.71) 
(9.72) 

(9.73) 
(9.74) 

(9.75) 



This means that the excitations of (j) cannot appear in physical states, while an arbitrary excitations 
of (all having zero norm) can accompany physical photons. The Hamiltonian which effectively 
acts on physical states has the form 



Lor 



phys 



(Tx 



1 ^ y 



+ 



16LL 
1 e2 



^(P-)^ + / d'xCl^A^) [ ^Q__\yphys * W) 



(x) 



+ 2 2L ' ^ x_i_rihy,^{x±) [g(±)[Mphys] * ^Ihys) ix±) 



(9.76) 



and is identical with the Hamiltonian ( |8.84| ) obtained for the LF Weyl gauge after the implementation 
of Gauss' Law. Thus we can take the above gauge-independent expression as the physical Hamiltonian 
for QED at the finite volume LF. 
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10 Transverse Coulomb gauge 



When the Coulomb condition is imposed for transverse components of gauge potentials ^jV^ = it 
effectively introduces constraints for two sectors 

Therefore we have to add some gauge condition for the global zero modes and we choose the LF Weyl 
gauge g+ = 0. Because all these conditions are non-dynamical, we will impose them explicitly via the 
suitable decomposition of vector gauge fields into transverse and longitudinal parts. 

10.1 Proper zero mode sector 

As usually, we start with the Lagrangian density 

^cZl = dia^dia+ + ^ (d+a^f - ^ [diaJ^'Y + J^zm(^+ + Jl^mf^J'' + Jpzm^^- , (10.2) 
where the transverse components of fields are defined as 

ar^(6ik-didk^*)ak . (10.3) 



The Euler-Lagrange equations 



'A 



dla- = -Aa+ + J-^ , (10.4) 
-did+a- = Axaf + 4,^, (10.5) 
= -Axa_ + J+„, (10.6) 

can be explicitly solved in terms of external currents 

«- = -^*Jpzm, (10-7) 

ai = * Jpzm + * * djJ^,^^ , (10.8) 

a+ = Jpzm-9+^* (^-^* J^,m^ . (10.9) 

So here again there are no independent proper zero modes and there is no canonical structure in this 
sector. Thus we only write the effective Lagrangian 

^pzm eff ^ j+ J_ J- _}_ ^ n _}_( J_ ^ a. ji 

Coul '^pzm '^pzm 2 '^pzm 2 ^I'^pzm 

- ^{Apzmf-d+Apzm^*J+^, (10.10) 

where the subsidiary non-dynamical field Ap^m has been introduced to avoid the explicit presence of 
(9, 7+ 

^+"pzm' 



64 



10.2 Normal mode sector 

Normal modes are described by the Lagrangian density 



CZni = diA+d,A^+d.Ard+Ar + ^id+A^-d.A+)' 



y + A^r + A+j+ + AJ^jZ" , (10.11) 



2 

which generates the Euler-Lagrange equations 

(2a+a_ - A±) A. = 9_ + d-A+) - J+„ , (10.12) 

(2a+a_-A^)A+ = d+ {d+A^ + d-A+) - J-^ , (10.13) 
(2a+5_-Aj_)yiP = Jjr. (10.14) 

These equations describe dynamical evolution but contain also the constraints. To see this most easily 
we parameterize two gauge field potentials A± by the single field $ 

A_ = , (10.15) 

A+ = (^A^cD-^* j+^j +5+$ , (10.16) 

which satisfies the dynamical equation of motion 

(29+9_-A^)$ = + J_*(5_J~^ + a+J+J • (10.17) 

In the canonical analysis the last equations would need a modification for removing the term (9+J+„. 
However this would introduce such a term in the definition of A^, therefore, just like in the case of 
massive QED in Section |l], there is no unconstrained canonical structure consistent with the primary 
constrained system. Thus all we can find is the effective Lagrangian density which contains dynamical 
fields ^ and Aj'' 

^Tout" = - \ (^A^a> - ^ * + 9_$J-^ + 9+$J+^ 

+ d^Ard^Aj^-\(d,Af)\Aj^jl::. (10.18) 



Combining two effective Lagrangians (10.10) and (10.18) we obtain the total effective Lagrangian for 
the complete system with fermion fields 

C^i, = d+A^^d.A^^ -\[d,A]^)\d+d,^d^d,^-\[A^^-^*Jt, 

- MV'+t7V--MV'-VV'+-e\/2^+V+^+°"'-e\/2V'_V-"^^-^°"' (10.19) 



2 
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where 

yCoul ^ ^Tr ^ ^. ^ (^0 20) 

yCoul ^ 5_^ + ^___L,j+^ ^ (10.21) 
yCoul ^ _ Q^±_ ^ ^^^^ _ (1Q_22) 

10.3 Fermion field sector 

Before writing down the Lagrangian for the fermion field we have decided to introduce the dressed 
fermion field 

i'±{x) = '4^±{x) exp * Ap2m(x_L)| (10.23) 

and this leads to 



— m 



where 

y| = d+^. (10.25) 

Now we have the same form of constraints for dependent fermions as in the case of LF Weyl gauge 
in Section ^ 

(iV2d. - eV2y5°"') ^- = -ia'di4}+ + Mj°i;+ + ea' (vf"''^ - ^ * -0+ , (10.26) 

= y d^x {^la'ip- +i^la'^+){x) , (10.27) 

thus we can adopt those results directly and only the equation for Ap^m is new 

A±Ap,^ = -dj;,^ . (10.28) 

We easily obtain 

■^4.m - % = -^S±[^'] * (10-29) 
and further the effective Lagrangian density for fermions takes the following form: 

(10.30) 

where 

C(x) = \Mj^-ia'di + ea'Vf'"'\x)\ip+{x), (10.31) 
f'(^^) = ^ / dx-^^(x)a- (^ ^^_ * (x) 
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1 



1 



id- 



The equation for the subsidiary dependent field A. 



pzm 



A, A, 



pzm 



(10.32) 
(10.33) 

(10.34) 



can be uniquely solved but we really do not need this solution. Since Apzm has disappeared from the 
effective Lagrangian ( |10.3C| ), we can omit it in further steps. 

Our analysis of constraints for the complete QED in the Coulomb gauge leads to the final effective 
Lagrangian density 



rCoul 



1 



1 e2 
2 



(10.35) 



10.4 Canonical quantization 



From the classical effective Lagrangian ( 10.35 ) we can infer the structure of quantum (anti) commutation 
relations for all independent fields. Besides the canonical relations 



2a! 



-i6,,Aj^6l^\x-y) , 



Ar{x),A,Ar{y) = 

2dlMx),Hy)] = iA^6^^\x-y) , 
{^P'Ux),^\{y)} = AW5(3)(x-y) , 

we find also two noncanonical commutators between $ and the fermion fields 

29^ [$(X([,x),V^+(a;+,y)] = e5^(x-y), 
25^ [$(x+ x),V^^(x+,y)] = -eS^{x-y), 

which come from the term The quantum Hamiltonian is assumed in the form: 

2' 



(10.36) 
(10.37) 
(10.38) 
(10.39) 

(10.40) 
(10.41) 



h: 



Coul 



quan 



J d'x \ (diA]^)' + \ (a^<J. - e^/2^ * (V^tV^+))" 



+ 



16LL 
1 

2 2L 



(10.42) 
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tacitly defining proper ordering of non-commuting propagators. Thus we see that the transverse 
Coulomb gauge condition has the same number of independent excitations as the physical subsystems 
defined either in the LF Weyl or in the Lorentz covariant gauges. Here one encounters a strange 
phenomenon that one boson field does not commute with fermion fields and this is connected with 
the rather extravagant choice of gauge which imposes the constraint on transverse gauge fields which 
are known to describe two physical photon excitations at LF. This non-commutativity is the lowest 
price we can pay for such careless choice of gauge fixing condition and it can be solved if we introduce 
another dressed physical fermion fields 



phys 



expie$(af)V'+(x) 



(10.43) 
(10.44) 



which already commute with $ fields. Now the quantum Hamiltonian ( 10.42 ) takes the form 



Tjphys 
quart 



2 1 

+ 2 



1 \' 



where 



+ 



+ 



phys 



16LL 
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2 2L 



(x) 



I Ml 



^phys ix±) , 

AT'{x)-d^^x)+qi{x+) , 



^physix) 
^phys (-^-L ) 



[X 



1 



2L 



(10.45) 

(10.46) 
(10.47) 
(10.48) 



—= / dx v. 



+ 



1 

72 
1 
71 



1 



phys 



X ]a 



Mlhy six !_) 



dx (physi^)^^" 

72 / dx-^l^y^{x) { 



1 



^id- — eV: 
1 



* Cphys j {x) 
phys * '^Phys I (^) 



^phys 



* ^phys 



;i0.49) 



;io.5o) 



Jd--eV!l 

and agrees with the former expressions for the physical Hamiltonians in the LF Weyl gauge ( ^.84| ) 
and the Lorentz covariant gauge (9.76). 



10.5 Physical gauge conditions 

The physical field variables that we have just found after redefinitions of the fields can be obtained 
directly from the following gauge conditions: 



Q+ 



, 

d±a± 
. 



(10.51) 
(10.52) 
(10.53) 



Only the sector of normal modes has not been analysed by us so far |44], therefore below we sketch 
some crucial points. From the Lagrangian density 



(d+A, - diA+) d-Ai + ^ (d.A+f - - 



^ {diAj - djAif + A+J+^ + , (10.54) 
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where the gauge fixing condition = has been exphcitly implemented, we derive the Euler- 
Lagrange equations 

a_ {d-A+ - d^Ai) = J+„ , (10.55) 

i2d+d- - Ax) A^ = di id^A+ - djAj) - 4^ . (10.56) 



Evidently, A^ is not a dynamical variable and its equation of motion (10.55) can be solved as 

A+ = ^*{9^^^ + ^*Jnm)- (10.57) 

Only one Dirac bracket is nonzero 

2d^_{Ai,{x),Aj{y)}^^ = -Sij6^{x-y) (10.58) 

and the Dirac Hamiltonian is 

Hd = lidiAi + j- * J+J' + i {d,A, - d,A,f - AJL ■ (10.59) 

These results can be incorporated into the whole procedure of equivalent Lagrangians if the above 
results are transformed into the effective Lagrangian density for the normal mode part 

C:]} = d+Ad.A, - ^{d,Ai + j- * 4+^)2 - 1 {d,A, - d,A,f + A^Jl^ . (10.60) 



Then following the analysis from Subsection 9.2 and (with trivial modifications) from Subsection [97 
one finds the commutator relations for physical fields ( |10.36 - 10.391) and the physical Hamiltonian 
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Part V 

Conclusions and perspectives 

In this paper we have presented a novel method of canonical quantisation for constrained systems 
when applying it to the analysis of QED at LF. Separating different sectors of fields, we have been 
able to find effective description of them avoiding an explicit implementation of Dirac's procedure for 
the whole system of entangled constraints. Different gauge conditions proved to be legitimate choices 
for the quantum gauge field system interacting with fermions. 

Feynman rules for perturbative calculations have been found and they have common properties. Gen- 
erally, the canonical propagators for gauge fields have additional noncovariant terms which behave 
non-causally i.e. in the Fourier representation they are given solely by the CPV poles. These terms 
are cancelled by the direct interaction of currents which canonically appear in the interaction Hamil- 
tonians. The fermion parts of interaction Hamiltonians are factorized and the non-covariant factors 
cancel with the additional non-covariant term in the fermion field propagator at LF. In this manner 
there are two equivalent settings of Feynman perturbative rules, the first (canonical) rules can contain 
non-causal and non-covariant terms, the second (effective) rules are causal and covariant. 
The LC-gauge has been investigated as the limit of the general axial gauge and the flow covariant 
gauge. The ML-prescription for the LC-gauge propagator has been derived only in the second choice 
and these results hold for the interactions with fermions. 

In the DLCQ method three choices of gauge conditions have been canonically analysed. In the cases 
of the LF Weyl gauge, the Lorentz covariant gauge and the transverse Coulomb gauge, the same 
physical Hamiltonian has been found though via different methods of Gauss' law implementation, 
the Lautrup-Nakanishi condition and the redefinition of fermion fields, respectively. Finally, the set 
of physical gauge conditions, which straightforwardly leads to the above physical Hamiltonian, has 
been proposed. 

All the above results have been derived for fermion currents. The problems which can appear for 
other charged matter fields have been discussed in the simplest model of 1+1 dimensional LF-Weyl 
QED with scalar fields. A reasonable analysis of this model has been formulated starting with the 
nonlocal redefinition of one scalar field at LF. For these new fields, the canonical analysis and struc- 
ture of perturbative calculations have been studied. Both methods and results were similar to those 
for fermion fields. 

We have found that the presence of field derivatives in the matter currents (as in the scalar case) 
can be neutralized at LF by a suitable nonlocal redefinition of these fields. We expect that this is a 
general property of the LF formulation and can also be applied for more physically relevant models. 
Therefore wc plan to analyse, first, the QED for charged vector fields and then the non-Abelian in- 
teractions. Especially the latter case will ultimately settle the status of the LC-gauge at LF, first of 
all the question of proper prescription for spurious poles within the canonical LF formulation. 

Acknowledgements 

I thank Professor Dominik Rogula for his encouragement and many invaluable remarks during my 

research on this project and Dr Stanislaw D. Glazek for many discussions on different aspects of 
light-front physics. I am also indebted to so many people who have helped me to solve all problems 
that I have encountered while preparing this thesis. 

I express my gratitude towards Dr Alex Kalloniatis and Dr Rik Naus for their collaboration in our 
joint paper which was the starting point of my present research. 



70 



Appendices 



A LF Notation 
A.l Coordinates 

In 3+1 dimensions we define longitudinal coordinates = and take as the dynamical 

evolution parameter. We denote transverse components x± = by Latin indices (i, j, ...). 

Similarly we define components of any 4-vector and a scalar product of two 4- vectors decomposes as 

A ■ B = A~^B~ + A~ — A^B''. The metric has non-vanishing components = l,gij = —Sij. 

Partial derivatives are defined as d± = d/dx^, di = d/dx^. Tensor components are defined analo- 
gously e.g. T^^' = -^{T^^ ± T^^) and summation over repeated indices is understood. 
Also we introduce a vector notation for components x = {x~,x±) which lie on a light-front surface 
x+ = const, and for momenta associated with them k = {k-,k±). Scalar product of such 3-vectors 
decomposes as k ■ x = k^x~ — kiXi. 

In the D+1 dimensions the notation generalizes trivially with the number of transverse directions 
changing from 2 to d. 



A. 2 Dirac matrices 

The Dirac matrices 7^ satisfy anticommutation relation 

+ = '^9^'' , (A.l) 

where their components are defined analogously to coordinates e.g. 7^ = . Thus 7^ are 

nilpotent matrices (7 ) = 0. 
For the projection operators 

A± = 4=7°7^ = ^7^7^ , (A.2) 



V2' ' 2 



we have useful relations 



A±A± 


= A±, 


A+A_ 


= A_A+ 


A+ + A_ 


= 1 , 


7±A^ 


= A±7^ 


7±A± 


= A±7± 


7°A± 


= A^7° 


yA± 


= A±y. 



In many places we also use the standard notation 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 



7° = (A.IO) 
7V = a\ (A.ll) 
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B Green Functions 



B.l Feynman propagator functions 

In 1+1 dimensions we define noncovariant Feynman Green functions Ep{x) and Ep{x) 



df 



^ 



e(a;+)e- 



-ik-X 



e(-x+)e 



ik—x 



1 ( e(x+) ^ e(-x^ 



27r \ X — ie x + ie J 
df a;+ / e(x+) G(-x+) 



+ 



27r \ X — ie x + ze 
which can be also represented by the 2-dimensional Fourier integrals 



oo 



oo (27r 



2 k+ + ie sgn(A;_) ' 

-ik-x ^ 



(27r)2 + sgn(fc_)]2' 

In D+1 dimensions we define covariant Feynman Green function as: 



D%+\x) 



df 



= I 



°° d'^k_i 

-oo (2^'^ Jo 2tt 2k. 

dd+'ij^ ^-ik-{x-y) 



(27r)<^+2 A;2 + ie ' 

The Feynman propagator function for the massive fields in 3+1 dimensions are given by 



Af(x,m2) 



d^k 



-ik-x 



(27r)4 2A;+fc_-A;2 -M2 + ie ' 
thus the functions for massless fields can be defined as the limits 



Df{x) 
Ef{x) 



lim Ai^fx.M^) = i 



d^k 



„—ik-x 



(27r)4 2k+k-. -kl + ie' 



lim 



d 



dm? 



Ap(x, m 



d^k 



-ik-x 



(27r)4 {2k+k- - kl + ief 



In the main text we use the following property 



2dlDp\x) + EUxL)S''{x±) 



- ^"'V-oo(2^io 

/ + 

+ Q{-x+)[e^''-'' 



dk- 



(B.l) 
(B.2) 

(B.3) 
(B.4) 



(B.5) 
(B.6) 

(B.7) 
(B.8) 



27r 

^-ik-{x-y) 



1 e- 



Jo 



(27r)<^+2 A;2 + ie A;+ + ie sgn(A;_) ' 
and a similar property holds also for the massive case 

d^k e-''=-(^-2') 1 



2d'_AF{x,m') + E},{xl)6\x±) = - (Aj 



m 



(27r)^ k"^ — m? + ie k+ + ie sgn(A;_) ' 



(B.9) 



(B.IO) 
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The noncovariant Feynman propagator functions in 3+1 dimensions are given by 

/d^k g-ik-x 
(27r)4 il + a)k+k.-akl+ie ' ^^'^^^ 

(27r)4 [(l + a)fc+fc_-afei+i€]2- ^^"^^^ 

B.2 Integral operators in x~ direction 

We use the basic inversion of d- defined by the Fourier integral 

dkexp —ik-{x^ —y~) CPV— = -sgn(a;~ —y~)- (B.13) 

However, for the fermion and scalar fields it is more convenient to use another Green function 
(idJ)~^{x~ — y~) = —i{d^)~^{x~ — y~) = — zisgn(a;~ — y~)-, which under the complex conjuga- 
tion behaves like a Hermitian matrix 

\id.)-\x- -y-)\ = {id.)-\y- -X-). (B.14) 

Therefore, the Green function for the covariant derivative {id- — eAJ)~^ can be given in the terms 
of integral operators 

{id- - eA_)-^[x- ,y-] = J dz-{id-)-^{x- - z-)W-i[z- ,y-;a\ (B.15) 

= J dz-W-i[x-,z-;a^]{id-)-\z- -y-) , (B.16) 

where 

oo 

W-i[x-,y-;a] = .^(x" - y") + ^(a)^[a;-, y"], (B.17) 

k=l 

a[x-,y-] = eA-{x-){id-)-\x- -y~), (B.18) 
{ar[x-,y-] = ldz-{ar-'[x-,z-]^z-,y~]= ldz-a[x-,z-]{ar-'[z-,y-], (B.19) 

at[x-,y-] = e{id-)-\x- -y')A^{y-), (B.20) 
{d^r[x-,y-] = j dz-{a^r-\x- ,z-p[z- ,y-] = j dz'a^x' , z-]{a}r-\z- ,y-]. 

(B.21) 

One can easily check that under the complex conjugation these expressions behave as follows: 

{{aT[x-,y-])* = {a^r[y-,x-] (B.22) 
((at)-[x-,y-])* = {ar[y-,x-] (B.23) 

{W-i[x-,y-;a]y = W-i[y- , x~ -,(1^] (B.24) 
(W_i[x-,y-;at])* = W-i[y-,x-;a] (B.25) 

(^{id- -eA-)-^[x-,y-]y = {id- - eA-)-^[y- ,x-] . (B.26) 

If all two-argument expressions are treated as generalized (infinite dimensional) matrices and their 
convolutions are denoted by asterisks then one can introduce the self-explanatory matrix notation 

{id- - eA-)-^ = {id-)-^ * W-i[a\ = W-i[a'] * {id-)~^ , (B.27) 
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where 



W-i[a] = l + ^(a)\ (B.28) 

k=l 

a = eA^{id-)-^ , (B.29) 

(a)" = {a)''-^*a = a*(a)''-\ (B.30) 

at = e{id-)~^A^ , (B.31) 

(at)" = (at)"-i*at = at*(at)"-i. (B.32) 

Also one can check the fohowing property 

id-W-i[a^* = id- + eA_W-i[a^]* , (B.33) 

so the integral operator W_i[at] effectively transforms the covariant derivative into the partial deriva- 
tive 

{id- - e^_)W_i[at]* = id-. (B.34) 
Another simple calculation shows that 

{id_ -eA-){id-r'^ *W-i[a] = W_i[a] - a * >V_i[a] = 1 (B.35) 

and this constitutes the proof that the integral operator defined by ( [B.27 ) is the Green function for 
the covariant derivative operator {id- — eA-). All above results can be easily transformed to the 
quantum theory, where the real- valued field A-{x^) is substituted by the Hermitian operator and 
the complex conjugation is replaced by the Hermitian conjugation. 

All charged fields can be easily incorporated into the matrix notation, where and the scalar field 4> and 
fermion field ^ are treated as the one-column matrices while their Hermitian conjugated counterparts 
^t and as the one-row matrices, respectively. 



B.3 Inverse Laplace operators 

The Green function for the Laplace operator is well defined for d > 2 dimensions 



and for (i = 2 it is singular. From we know that there are also other regularizations, strictly in 
d = 2 dimensions. In the paper, we use the massive regularization - when the pole aX k\ = Q is shifted 
by the mass parameter 

[AJ-i(xJ - [A^-mr\x^) = -J ^^Jj^ ^ = --Ko(mV4), (B.37) 

which naturally appears for the massive electrodynamics. 
Another regularized Green function appears in Section |5| 

1 r d^ki dk- 



e 



-ik-x 



and for a < there is no ambiguity in the integrand. 
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B.4 Finite volume Green functions 

One can safely invert differential operators such as (9„ and d'j^ in terms of well-defined Green's func- 
tions, taking care of the respective mode sector in which the operator acts. With the covariant 
derivative = 9^ + ieV^, we define the operator- valued Green's function to [iD-): 

(iD^. ) .„ ^ [x, y\ ^ 5^'^ (x-y)- [Sub.] , (B.39) 
lO- — eH 

where [Sub.] denotes possible subtractions corresponding to zero eigenvalues of the operator in ques- 
tion. A nonperturbative construction for this particular Green's function is presented in Appendix A 
of [40|. The Green's function G{±)[x±,y±; O] is defined by the relation 



e2 



[Ai - ^ 6^^\x^-y^) , (B.40) 



where A_|_ = d\ , and now O is some field operator of mass dimension one. Eg. ( B.40 ) can be elucidated 
order by order in perturbation theory for which one uses the basic inversion of A_l. A nonperturbative 
definition is however nontrivial and therefore the above operation is, at best, merely formal and its 
concrete implementation remains an open problem. 

The convolutions of the above Green's functions are also denoted by asterisks. 
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